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ABSTRACT 

We present an optical spectroscopic study of a 90% complete sample of nearby ULIRGs 
(z < 0.175) with optical Seyfert nuclei, with the aim of investigating the nature of the 
near-nuclear (r < 3.5 kpc) warm gas outflows. A high proportion (94%) of our sample 

show disturbed emission line kinematics in the form of broad (FWHM > 500 km s _1 ) 
and/or strongly blueshifted (AV < —150 km s _1 ) emission line components. This 
proportion is significantly higher than found in a comparison sample of nearby ULIRGs 
that lack optical Seyfert nuclei (19%). We also find evidence that the emission line 
kinematics of the Sy-ULIRGs are more highly disturbed than those of samples of non- 
ULIRG Seyferts and PG quasars in the sense that, on average, their [OII]AA5007,4959 
emission lines are broader and more asymmetric. 

The Sy-ULIRG sample encompasses a wide diversity of emission line profiles. In 
most individual objects we are able to fit the profiles of all the emission lines of different 
ionization with a kinematic model derived from the strong [OIII]AA4959,5007 lines, 
using between 2 and 5 Gaussian components. From these fits we derive diagnostic 
line ratios that are used to investigate the ionization mechanisms for the different 
kinematic components. We show that, in general, the line ratios are consistent with 
gas of super-solar abundance photoionized by a combination of AGN and starburst 
activity, with an increasing contribution from the AGN with increasing FWHM of 
the individual kinematic components, and the AGN contribution dominating for the 
broadest components. However, shock ionization cannot be ruled out in some cases. 
Our derived upper limits on the mass outflows rates and kinetic powers of the emission 
line outflows show that they can be as energetically significant as the neutral and 
molecular outflows in ULIRGs - consistent with the requirements of the hydrodynamic 
simulations that include AGN feedback. However, the uncertainties are large, and more 
accurate estimates of the radii, densities and reddening of the outflows are required 
to put these results on a firmer footing. 
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1 INTRODUCTION 

AGN-induced outflows are increasingly recognized as a key 
element in the overall galaxy formation process, regulating 
the correla tions between black hole mass a nd host galaxy 
properties (|Silk fc Reesl . Il998l : iFabianl . 119991 ). and influenc- 
ing the luminosity fun ctions of galax i es, pa rticularly at the 
high luminosity end (|Benson et al.l . [2003). Indeed, AGN 
feedback is now routinely incorporated into numerical sim- 
ulations of the hierarchical growth of galaxies through ma- 

* E-mail:javierr@iac.es 



jor g alaxy mergers (|Springel et all . 120051 : | Johansson et alj , 
120091 ). However, there remain considerable uncertainties 
about the exact nature of the symbiosis between black holes 
and their host galaxies. Consequently, the importance of 
AGN-induced outflows relative to those driven by the su- 
pernovae associated with the merger-induced starbursts has 
yet to be established observationally. 

In order to test the models it is important to identify 
samples of actively evolving galaxies in which both the black 
holes and the galaxy bulges are growing rapidly. Although 
much attention has been paid recently to sub-mm galax- 
ies at high redshifts that frequently show signs of AGN ac- 
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tivity in X-ray obser vations ([Alexander et all . 120051 . |2010| ; 
lHarrison et all [20121 ). there is a limit to what can be learnt 
in detail about the co-evolution of black holes and galaxy 
bulges in such objects because of their faintness. Moreover, 
for high redshift galaxies many important diagnostic emis- 
sion lines are shifted out of the optical/near-IR, and their 
outflows are difficult to resolve spatially. 

Fortunately, the hierarchical evolution of galaxies is 
also continuing in the local Universe, albeit at a reduced 
rate. The Ultra Luminous Infrared Galaxies (ULIRGs: 
Lir > 10 12 L0) represent a class of objects at rela- 
tively low redshifts in which the prodigious far-IR radia- 
tion represents the dust re-processed light of major star- 
bur sts and /or AGN buried in the nuclei of the galax- 
ies ((Sanders fc Mirabell . 1 19961 ). ULIRGs almost invariably 
show morphological evidence (e.g. tidal tails, double nuclei) 



mergers (e.e. ISanders et a 


11 19881: ISurace et all 1 19981. 200C: 


Surace & Sanders! 12000; 


Kim et all 20021; Veilleux et al. . 


2002). Based on their mid- to far-IR colours they can be 



sub-divided into "warm" (/as^m/ZeOfim > 0.20) and "cool" 
(/25/jm//60nm < 0.20) sub-types, with the warm sub-type 
freq uently associated with v isible Seyfert or quasar nuclei 
(e.g. ISurace fc Sanders! [1999). In terms of evolutionary sce- 
narios, it has been proposed that cool, starburst-dominated 
ULIRGs evolve into warm, AGN-dominated ULIRGs as the 
circum-nuclear d ust is dispersed by starburst and AGN - 
induced outflows (jSanders et all 1 19881 ; ISurace et all ll998T l . 
Such scenarios are broadly consistent with current simu- 
lation s of major, gas-rich galaxy mergers l|Hopk ins et all 
2005). Therefore, the ULIRGs represent local analogues of 
the rapidly evolving galaxies now being detected in large 
numbers in the distant Universe, but they have the consid- 
erable advantage over the high redshift objects that they 
are close enough to study in detail. In this sense the nearby 
ULIRGs are key objects for understanding the detailed phys- 
ical mechanisms - particularly the outflows - involved in 
the co-evolution of super-massive black holes and their host 
galaxies. 

As the most rapidly evolving galaxies in the local Uni- 
verse, ULIRGs represent just the situation modelled in many 
of the most recent hydrodynamic simulations of gas- rich 
mergers (jdi Matteo et al ], 120051 : IJohansson etHl l2009h . In 
the final stages of such mergers, as the nuclei coalesce, the 
simulations predict that the super-massive black holes ac- 
crete at sufficiently high rates to produce luminous, quasar- 
like AGN; the outflows associated with the AGN are poten- 
tially capable of ejecting the remaining gas and halting star 
formation in the bulges of the merger remnants. However, 
fine tuning of the models is required in order for them to 
reproduce the observed M^h vs. a relationship. In particular 
it is required that the outflows carry a relatively large pro- 
portio n of the avai l able accretion power o f the AGN (~5 - 
10%: iFabianl . Il999l ; Idi Matteo et~all 12005! but see Hopkins 
& Elvis 2010). Clearly it is important to use direct observa- 
tions to determine whether this fine tuning is justified: are 
the AGN-driven outflows in ULIRGs truly as important as 
many of the simulations require? 

In line with the expectations of the hydrodynamic 
simulations, spectroscopic evidence for outflows has now 
been found in several ULIRGs usi ng optical ob s ervatio ns 
of neutral NalD absorption lines (|Rupke et all [2005a b; 



iRupke fc Veilleux! I201l1 ). and far-IR/sub -mm observations 
of the molecular OH and CO lines jFeruglio et all |201()I ; 
iFischer et all l20ld : ISturm et all l201lh . However, the ev- 
idence for a clear link between these outflows and nuclear 
AGN activity is mixed. For example, although the mass out- 
flow rates and kinetic powers associated with the neutral 
NalD absorption line outflows are significant, comparisons 
between the properties of the neutral outflows in objects 
with an d without opt i cal AG N show surprisingly little dif- 
fe rence jRupke et al. . [2005 J ). Moreover, while the results 
of ISturm et al.l ( 2011 ). based on Herschel data, do suggest 
that the molecular outflows may be fas ter and more power- 
ful in ULIRGs hosting powerful AGN (|Sturm et all hoill ). 
the sample is small (only 6 objects) and incomplete. 

It is important to recognize that the gas associated with 
any AGN-driven outflow may be highly ionized and therefore 
not adequately sampled by NalD absorption and molecular 
line studies. In addition, the quantitative estimates of the 
properties of the neutral and molecular outflows are com- 
promised by the fact that, in most cases, the radial extents 
of these outflows are not well constrained by the observa- 
tions, which sample gas along the line of sight. 

Results based on the optical studies o f the large-scale 



(>5 k pc) warm ionized ga s in U L IRGs ([Heckman et al 
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reveal a similar picture to the molecular 
and neutral gas studies: energetically significant outflows, 
but no clear differences between the outflows in ULIRGs 
with and without AGN; again it is not clear whether the 
large-scale emission line gas is truly sampling the AGN- 
driven outflows. However, optical and mid-IR spectroscopic 
observations have recently revealed highly ionized out- 
flows in the near-nuclear regions of several ULIRGs with 
Seyfert-like nuclei in the form of b road, blueshifted emission 
line components (IWilman et all 1 1999; Li'p ari et all 120031 ; 
[ Holt et all 120031 . 120061 ; ISpoon et all 12009 ; ISpoon fc Holtl 
12009! ) . These observations also sugges t that the near- 
nuclear outflows ar e highly stratified (|Holt et al.! 120031 ; 
ISpoon fc Holtl . 120091 ). and that relativist ic jets may play 
an important r ole in accelerating the gas (|Holt et all 120091 ; 
ISpoon fc Holt] . 120091 ). Given that they isolate the AGN- 
driven outflow components, the near-nuclear high ionization 
emission lines thus provide a powerful tool for addressing the 
key questions surrounding AGN-driven outflows in ULIRGs. 
These questions include the following. 

• What is the dominant outflow driving mechanism in 
ULIRGS: AGN or starbursts? As discussed above, previous 
studies based on observations of the molecular, neutral and 
large-scale ionized gas outflows have not provided a clear-cut 
answer to this question. Clearly it is important to consider 
this question in the context of the near-nuclear warm out- 
flows. 

• Are AGN-driven outflows more significant in Seyfert- 
like ULIRGs than in other types of AGN? Given that the 
ULIRGs are among the most actively evolving galaxies in 
the local Universe, with (predicted) high nuclear nuclear ac- 
cretion rates, we might expect the AGN-induced outflows 
in such objects to be unusually strong. However, a large 
proportion of AGN of all types show evidence for outflows 
at some level, in the form of blue wings to the [QUI] emis- 
sion lines fe.g. lHeckman et all Il98ll : IWhittleL 1 19851 ; IVeilleuxl . 
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Il991bl ; iNeTson fc Whittle . Il995l ). Therefore, it is important 
to establish whether the outflows in Seyfert-like ULIRGs are 
significant more extreme than those detected in the general 
population of non-ULIRG AGN. 

• What is the nature of the near-nuclear outflows in 
ULIRGs? Comparisons of spectroscopic observations of 
emission lines of different ionisation have provided evi- 
dence that the near-nuclear emission lines outflows in some 
ULIRGs are stratified. However, it is currently uncertain 
whether the outflows comprise a small number of distinct 
outflow velocity components, each associated wit h a particu- 
lar sp atial location and set of physical condition l)Holt et all 
2003), or rather encompass a continuous gradient of veloc- 
ities and physical condit ions across the full extent of th e 
emission line region (e.g. iDe Robertis fc OsterbrockL Il986h . 
At the same time it is important to confirm using a wide 
range of emission lines that the outflows are indeed deceler- 
ating, as sugges ted by observat i ons o f the mid-IR fine struc- 
ture neon lines (|Spoon fc Hohl 120091 ) . 

• Do relativistic jets play an important role in accelerat- 
ing the ULIRG outflows? Although it is certainly the case 
that some of the most extremely blueshifted emission line 
components are are found in ULIRGs t hat contain powerful 
relativistic jets (jrlolt et all 120031 . l2006n . the issue of the jet 
contribution to the acceleration and ionisation of the out- 
flows warrants further investigation using a complete sample 
of ULIRGs that is unbiased towards jet properties. 

• How energetically significant are the AGN-induced out- 
flows in ULIRGs? In terms of testing the hydrodynamical 
merger simulations that incorporate AGN feedback as a key 
element, it is crucial to accurately quantify the mass outflow 
rates and kinetic powers of the near-nuclear outflows. In this 
way it will also be possible to compare the properties of the 
near-nuclear AGN outflows with those of the neutral and 
large-scale ionised gas outflows, to determine whether the 
AGN-induced outflows are truly energetically significant. 

In this paper we directly address the these questions us- 
ing detailed optical spectroscopic observations of a complete 
sample of 17 nearby (z < 0.175) ULIRGs with Seyfert-like 
nuclear spectra. Throughout this paper we assume a cosmol- 
ogy with H = 71 km s _1 Mpc -1 , Q = 0.27, Ov = 0.73. 



2 SAMPLE SELECTION, OBSERVATIONS 
AND DATA REDUCTION 

In order to achieve the aims of this project we require a 
sample of ULIRGs with Seyfert-like nuclear spectra that are 
close enough to study in detail and have high quality nuclear 
spectra; the sample also needs to be complete as possible. 
As a starting point we t o ok th e "Extended Sample" (ES) of 
IRodnguez Zaurfn et al.l (|2009l ). which comprises 36 nearby 
[z < 0.175) ULIRGs o f all sp ectral types from the Uy sam- 
ple of iKim fc Sanders! i] 19981 ) that have deep WHT spectra. 
We then selected all the objects with Seyfert-like nuclear 
spectra fr om this sample. T o make this selection we use the 
criteria of lYuan et al.l (|2010T ). which in turn are ba s ed on the 
emission line diagnostic diagrams o f lKewlev et ail (120061 ) for 
SDSS galaxies. Specifically, we selecte d all the objects clas - 
sified as Sy2 in at least 2 out of 3 of the lKewlev et al.l (|2006l ) 
diagnostic diagrams, as well as the three objects with broad- 
line Seyfert 1 spectra. This results in a sample of 19 nearby 



ULIRGs with Seyfert-like spectra, redshifts z < 0.175, decli- 
nations 5 > —25, and right ascensions 12 < RA < 02 hr (see 
Table 1 for details). In co mparison, a tota l of 22 ULIRGs 
from the full Uy sample of IKim fc Sanders! (Il998l ) meet the 
same selection criteria in optical classification, redshift, dec- 
lination and right ascension, and of these only 3 objects - 
F12265+0219 (3C273), F13443+08, F14121-0126 - are not 
included in our sample of ULIRGs with Seyfert-like nuclear 
spectra and deep WHT spectra. Therefore our sample is 
86% complete^, 

and is certainly representative of the subset 
of ULIRGs with optic al Seyfert nuclei in the l o cal Un iverse. 

As described in IRodnguez Zaurfn et al.l (120091 ). deep 
spectra of all of the 19 objects in our sample were taken 
with the ISIS dual-beam spectrograph on the 4.2-m William 
Herschel Telescope (WHT), on La Palma, Spain. The obser- 
vations were carried out using the R300B grating with the 
EEV12 CCD, and the R316R grating with the MARCONI2 
CCD on the blue and the red arm respectively. The instru- 
mental setup resulted in a spatial scale of 0.4 arcsec/pix for 
both arms, and a dispersion of 1.72 A/pix for the blue and 
1.65 A/pix for the red arm. The useful wavelength range 
is ~ 3300 - 7800 A, in the observed frame. A slit of width 
1.5 arcsec was used for all the objects. In order to minimize 
the effects of differential atmospheric refraction, the objects 
were observed either with the slit aligned along the paral- 
lactic angle, or at small airmass (AM < 1.1). Although we 
lack accurate seeing estimates for individual objects, because 
ISIS/ WHT does not provide imaging observations, measure- 
ments by the DIMM seeing monitor at the observatory sug- 
gest that the seeing was in the range 0.6 < FWHM < 2.6 
arcsec, and usually less than 2.0 arcsec, for the nights of the 
observations. 

The data were reduced (bias subtracted, flat field cor- 
rected, cleaned of cosmic rays, wavelength calibrated and 
flux calibrated) and straightened before extraction of the in- 
dividual spectra using the standard packages in IRAF and 
the STARLINK packages FIGARO and DIPSO. The wave- 
length calibration accuracy, measur ed as the mean shift be - 
tween the measured and published (jOsterbrock et all . Il996l ) 
wavelength of night-sky emission lines, is ~ 0.35 A for the 
blue spectra and ~ 0.25 for the red spectra. The spectral res- 
olutions, calculated using the widths of the night-sky emis- 
sion lines (FWHM), are in the ranges 4.9 - 5.6 A and 5.1 
- 5.5 A for the blue and red arms, respectively. Finally, 
the estimated uncertainty for the relative flux calibration is 
±5%, based on comparison of the response curves of several 
spectrophotometric standard stars observed during all the 
observing runs. 

Following reduction, inspection of the spectra revealed 
that in 2 of the 3 objects with Seyfert 1 spectra (Mrk231, 
F21219-1775), the strength the broad Balmer and Fell emis- 
sion features and non-stellar continuum prevented analy- 
sis of their narrow line emission. Therefore we excluded 
these two objects from further analysis. Note however that, 
although it has not proved possible to study the narrow 
emission line kinematics in the nucleus of the two excluded 



1 This rises to 90% if we exclude from consideration the pow- 
erful quasar F12265+0219 (3C273), which is dominated by non- 
thermal jet emission, rather than thermal dust emission, at in- 
frared wavelengths. 
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Seyfert 1 objects, one of these objects - Mrk231 - is well 
known to show evidence for outflows based on absorption 
line studies of the neutral and molecular gas, and the ex - 
tended emission line regions (e.g. iRupke fc Veilleuxl . |201 ll) . 
Overall our final sample comprises 17 Seyfert-like ULIRGs 
with viable spectra. 



3 RESULTS 

In this paper we are interested in the AGN-induced outflows 
in the near-nuclear regions of ULIRGs. Therefore we concen- 
trat e on the 5 kpc-diameter nuclear ape rtures described in 
the iRodrfguez Zaurin et alJ l)2009l . |2010T ) stellar population 
study of the galaxies. In this section we describe the results 
obtained by fitting the profiles of the emission lines detected 
in the nuclear spectra. 

3.1 Fitting the emission line profiles 

Prior to the modelling of the emission line profiles in the 
nuclear regions, the spectra were shifted to the galaxy rest 
frame and a continuum model was subtracted from the data. 
In each case the continuum model was selected among the 
best fitting mo dels to the stellar continuum emission from 
the galaxy (see IRodrfguez Zaurin et all 120091 . for details on 
the stellar continuum modelling technique). 

Once the continuum emission was subtracted, we used 
the DIPSO package to fit Gaussian profiles to the emission 
lines. Our general modelling philosophy is to fit each emis- 
sion line with the minimum number of Gaussian components 
required to produce an acceptable fi10. The velocity widths 
derived from the fits were quadratically corrected for the 
instrumental profile, and all linewidths and radial velocity 
shifts were corrected to the rest frame of the objects. 

One of key aspects we want to investigate is the nature 
of the outflows: whether they comprise a small number of 
distinct outflow velocity components, each associated with a 
particular spatial lo c ation and set of physica l conditions (e.g. 
IVillar-Martfn et al.l . fl999l ; iHolt et "ail . l2008l l . or on the con- 
trary, there is a continuous gradient of velocities and physi- 
cal condi tions across the full extent of the emis sion line re- 
gion (e.g. lHolt et~alll2003l ; ls"poon fc Holtl.l2009D . In the first 
scenario one kinematic model would provide an adequate fit 
to all the different emission lines, while in the second sce- 
nario each line would have different kinematic components. 

To test these scenarios it is first necessary to select a 
prominent emission line that is suitable to generate a model 
that will be used to fit all the other emission lines. In this 
context the strong [OIII]AA4959,5007 lines usually have high 
S/N in our spectra, are in a region clean from any atmo- 
spheric absorption, and are not blended with other emis- 
sion lines (as in the case of Ha and [Nil] AA6548, 6583 or the 
[OII]AA3727,3729 doublet). Therefore, where possible, we 
started the individual analysis of each source by fitting the 
[OIII] AA4959,5007 emission lines. When fitting these lines 
we used three constraints in accordance with atomic physics: 

2 Wc define a "acceptable" fit to a emission line based on vi- 
sual inspection, i.e. one that has acceptable low residuals and 
adequately represents the overall shape of the line profiles. 



i) the flux ratio between [OIII]A4959 and [OIII]A5007 was 
set at 2.99:1 (based on the transition probabilities); ii) the 
widths of the corresponding kinematic components of each 
line were forced to be equal; iii) the shifts between the cor- 
responding components of each line were fixed to be 48. OA. 
We will refer to this as 'the [OIII] model' hereafter. 

We then attempted to model the other prominent emis- 
sion lines in the spectra with the same kinematic model (ve- 
locity widths and shifts) as [OIII] , leaving the relative fluxes 
in the kinematic sub-components to vary. As well as the con- 
straints derived from the [OIII] model, it is possible to fur- 
ther constrain the fits to other doublets in accordance with 
atomic physics: the shifts between the corresponding com- 
ponents of each line in all doublets (e.g. [NeIII]AA3868,3968, 
[NeV]AA3346,3425, [NII]AA6548,6583, [OI]AA6300,6363 and 
[SII]AA6716,6731) were set, and for some doublets (e.g. 
[Ne III]AA3868,3968, [Ne V]AA3346,3425, [NII]AA6548,6583, 
and [OI] AA6300,6363) the flux ratios were set based 
on the transition probabilities. For [SII] AA6716,6731, the 
flux ratio was required to be within the range 0.44 < 
[SII] (6716/6731) < 1.42, the ratios corresponding to the 
high and low density limits respectively. For simplicity, The 
[Oil] AA3726, 3728 lines were treated as a single line for the 
majority of the objects in our samplc[f]. Therefore, in those 
cases we assumed 3727A as the rest frame wavelength. 

The [OIII] modelling technique was successful for all 
but 4 of the 16 objects in the sample discussed in this 
paper (75%). This is consistent with the idea that, in 
the majority of cases, there exist a discrete number kine- 
matic components, each with its own line width, veloc- 
ity shift, physical conditions and ionization. However, in 
the cases of F13451+1232, F14394+5332, F16156+0146 and 
F23389+0300 (25% of the objects in our sample) it was not 
possible to use a single model to fit all the emission lines 
observed in their optical spectra, suggesting that there may 
be a continuous gradient of density, ionization and kinemat- 
ics through the Narrow Line Region (NLR), leading to each 
line having a different profile. 

Figure Q] shows the fits to the [OIII] emission lines pro- 
file for all but 3 of the objects in our sample. In the cases 
of F00188-0856 and F23327+2913, the [OIII]AA4959,5007 
emission lines have a low equivalent width and the profiles 
are affected by residual structure in the underlying contin- 
uum; hence it is not possible to precisely fit these lines. 
However, it is still possible in these cases to find one kine- 
matic model that fits all emission lines, based on the stronger 
Ha+[NH] lines. Figure[T]shows the fits to the Ha+[NII] com- 
plex for these two objects. In addition, a detail ed study of 
the ki nem atics for F1 3 451+ 1232 is presented in lHolt et al.l 
120031 ) and lHolt et~aH |201ll '). including the fit to the [OIII] 
emission line. Therefore, no fit for this source is shown in 
Figure[T] Table [2] presents the widths and shifts of the differ- 
ent kinematic components for the [OIII] model (or Hq+[NII] 
model in two cases) for the central 5kpc aperture for all 



3 There are 3 cases (F12072-0444, F17044+6720 and 
F17179+5444) for which accounting for the doublet nature 
of at least one of the different kinematic components, was 
necessary to adequately model the [Oil] AA3726, 3729 emission 
lines. 
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IRAS z RA DEC L m L[oin] Ll.4GH* 

FSC (J2000.0) (J2000.0) (W) (W) (W Hz" 1 ) 

(1) (2) (3) (4) (5) (6) (7) 



F00188- 


-0856 


0. 


.128 


00 


21 


26 


.5 


-08 


39 


26 


9 


.2E+38 


2 


1E+33 


6 


.8E+23 


F01004- 


-2237 


0. 


.117 


01 


02 


51 


.2 


-22 


21 


51 


7 


.5E+38 


1 


6E+34 


1 


.4E+23 


F12072- 


-0444 


0. 


.129 


12 


09 


15 


.4 


-05 


01 


14 


9 


.4E+38 


5 


.4E+34 


3 


3E+23 


F12112+0305 


0. 


.073 


12 


13 


16 


.0 


02 


48 


38 


8 


.2E+38 


1 


6E+33 


5 


.7E+23 


F13305- 


-1739 


0. 


.148 


13 


33 


15 


.2 


-17 


55 


01 


6 


.8E+38 


4. 


.4E+35 


2 


.6E+24 


F13428- 


f5608 


0. 


.037 


13 


11 


a 


.8 


55 


53 


14 


5 


.3E+38 


9 


.OE+33 


4. 


6E+23 


F13451- 


fl232E 


0. 


.122 


13 


47 


33 


.3 


12 


17 


24 


8 


.OE+38 


1 


2E+35 


1 


.9E+26 


F14394- 


f5332E 


0. 


.104 


14 


41 


04 


.3 


53 


20 


08 


1 


.7E+38 


1 


.1E+34 


1 


.OE+24 


F15130- 


-1958 


0. 


.109 


15 


15 


55 


.6 


-20 


09 


18 


5 


.3E+38 


3 


.7E+34 


2 


.7E+23 


F15462- 


-0450 


0. 


.100 


15 


48 


56 


.8 


-04 


59 


34 


6 


.2E+38 


3 


.4E+34 


3 


1E+23 


F16156- 


(-0146NW 


0. 


.132 


16 


18 


08 


.2 


01 


39 


21 


4 


.7E+38 


6 


OE+34 


3 


6E+23 


F17044- 


f6720 


0. 


.135 


17 


04 


28 


.4 


67 


16 


23 


5 


.8E+38 


1 


9E+34 


1 


.9E+23 


F17179+5444 


0. 


.147 


17 


18 


55 


.1 


54 


41 


50 


6 


.8E+38 


5 


.OE+34 


1 


.7E+25 


F23060- 


f0505 


0. 


.173 


23 


08 


34 


.2 


05 


21 


29 


1 


.2E+39 


1 


.OE+35 


5 


.OE+23 


F23233- 


f2817 


0. 


.114 


23 


25 


18 


.7 


28 


34 


19 


4 


.3E+38 


7 


.8E+34 


2 


.1E+23 
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f2913S 


0. 


.107 


23 


35 


11 


.9 


29 


30 


00 


1 


.9E+38 


5 


5E+33 


2 


1E+23 


F23389- 
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0. 
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23 


41 


31 


.1 


03 


17 


31 
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.3E+38 
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.7E+34 
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3E+25 



Table 1. The sample of 17 ULIRGs discussed in this papers. These objects were selected f rom t he iRodrfguez Zaurin et al. | Hooi) ES 
sample of 36 ULIRGs and are classified as Sy2 in at least two of the 3 of the lKewlev et all l|2006T l diagnostic diagrams. Col (1): object 
designation in the IRAS Faint Source Catalogue Database (FSC). For those sources with multiple nuclei, the individual nucleus considered 
for this study is indicated in the table (N = north, S = south, E = east and W = west). Col (2): optical redshifts from Kim &; Sanders 
(1998). Cols (3) and (4): right ascension (hours, minutes and seconds) and declination (degrees, arcminutes and arcseconds) of the IRAS 
source positions as listed in the Faint Source Database (FSDB). Col (5): IR luminosity from Kim and Sanders (1998) adapted to our 
cosmology. Col (6): total [OIII]A5007 emission line luminosities obtained from our modelling results (see text for details). Col (7): 1.4 
GHz monochromatic radio luminosities (Li.4GHz)- 



galaxies in our sample, whereas Table [3] presents the line 
fluxes and emission line ratios. 

We emphasize that the fits to some emission line blends 
have the potential to suffer from degeneracies, in the sense 
that it is possible for the fitting programme to allocate the 
fluxes between the different kinematic components in various 
ways and still produce acceptable fits to the overall emission 
line profile. This is a particular problem in cases where there 
are multiple, broad kinematic components, and the compo- 
nents of the blend are relatively close together in wavelength 
(e.g. Ha+[NII]AA6548,6584 and [SII]AA6717,6731). To check 
whether the degeneracy is a consequence of the particular 
fitting programme (DIPSO) we have used, we have also mod- 
elled the emission line profiles of those objects for which 
the degeneracies were more important ( e.g. F01004-2237 o r 
F15130-1958) using the MPFIT code |Markwardtl . l2009t) . 
The main difference between DIPSO and MPFIT is that, 
whereas DIPSO fits the kinematic components in each indi- 
vidual emission line or blend individually, we have used MP- 
FIT to fit all the kinematic components in all the lines/blend 
in the spectrum simultaneously. Overall, the modelling re- 
sults were consistent with those obtained using DIPSO and 
the degeneracies were still present. The impact of these po- 
tential degeneracies on the fits to emission line blends in 
individual objects is discussed in the Appendix. 

As well as fitting the 17 objects in our sample of 
ULIRGs with optical Seyfert nuclei, for the purposes 
of comparison, we used identical techniques to fit the 
[OIII]AA4959,5007 and/or Hq+[NII] profiles in the 5kpc 
aperture nuclear s pectra of the 16 nearby ULIRG s in the 
ES (z < 0.175) of iRodrfguez Zaurin et all l|2009h that do 
not show evidence for Seyfert nuclei at optical wavelengths. 



For ease of reference in the following sections, we use 
the following scheme to label kinematic components, based 
on line widths (FWHM): 

• narrow: FWHM < 500 km s" 1 ; 

• intermediate: 500 < FWHM < 1000 km s _1 ; 

• broad: 1000 < FWHM < 2000 km s" 1 ; 

• very broad 2000 km s _1 > FWHM. 

Detailed descriptions of the line profile fits for individ- 
ual ULIRGs are described in the Appendix, while in the fol- 
lowing sections we describe the general results derived from 
these fits. 
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Object FWHM AV 

IRAS 

km s _1 km s — 



(1) 


(2) 


(3) 


(4) 


F00188-0856 


N 


269±14 






I 


904±79 


-253±36 


F01004-2237 


N 


unres 






I 


849±79 


-229±32 




B 


1590±94 


-999±80 


F12072-0444 


N 


275±18 






I 


525±36 


-276±30 




B 


1343±84 


-446±35 


F12112+0305 


Nl 


162±73 


- 




N2 


473±66 


-171±14 


F13305-1739 


N 


435± 29 






Bl 


1275±31 


-36±14 




B2 


1685±286 


-281±29 


F13428+5608 


N 


450± 13 






B 


1368±66 


+14±11 


F13451+1232 


N 


340±23 






Bl 


1255±41 


-402±9 




B2 


1944±65 


-1980±36 


F14394+5332 


Nl 


335±26 






N2 


390±32 


-700±10 




N3 


435±52 


-1358±18 




Bl 


1272±67 


-156±8 




B2 


1401±131 


-1574±19 


F15130-1958 


11 


545±85 






12 


700±222 


-350±281 




B 


1630±42 


-725±131 


F15462-0450 


N 


147±57 






B 


1426±39 


-822±25 


F16156+0146 


N 


unres 






I 


804±26 


-186±10 




B 


1535±63 


-374±26 


F17044+6720 


N 


290±20 






B 


1765±103 


-553±65 


F17179+5444 a 


N 


358±75 


-123±33 




I 


515±33 


123±34 




B 


1562±43 


-242±61 


F23060+0505 


Nl 


376±10 


-152±6 




N2 


unres 


152±4 




B 


1001±23 


-463±20 




VB 


2150±125 


-1073±122 


F23233+2817 


N 


178±37 






I 


640±24 


-161±15 




B 


1511±118 


_447±40 




VB 


3433±903 


-2706±422 


F23327+2913 


N 


109±12 






I 


804±14 


-94±4 


F23389+0300 


N 


289±12 






VB 


2223±30 


47±13 



Table 2. The FWHM and velocity shifts (AV) for the different kinematic components of the galaxies in our sample. Col (1): object 
name. Col (2): The label of the different components as defined in the previous section [N (narrow), I (intermediate), B (broad) and VB 
(very broad)]. For those object with two components within the same FWHM range, these are indicated with numbers (e.g. Bl and B2 
corresponding to the two broad components detected in IRAS F13305-1739). Col (3) and (4): rest-frame widths (FWHM) and shifts 
relatively to galaxy rest frame (AV) of the different components. 

a Two models adequately reproduced the emission lines in the case of F17179+5444. One model comprises a narrow, an intermediate 
and a broad components while the second model includes two narrow components plus a broad one. The results show in the table for 
this galaxy correspond to the first model. For details on the modelling for this and all the other objects in our sample, see Section 3.2 
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IRAS F00188-0856* IRAS F01004-2237 




6500 6550 6600 6650 4900 4950 5000 5050 

Wavelength (A) Wavelength (A) 



IRAS F12072-0444 IRAS F121 12+0305 




'4900 49S0 5000 5050 '4900 4950 5000 5050 

Wavelength (A) Wavelength (A) 



IRAS F13305-1739 IRAS F13428+5608 
3.4 I 1 1 1 1 1 1 1 1 1 1 1 1 1 1 , 1 1 1 1 1 1 1 1 1 1 1 — 




'4900 4950 5000 5050 4900 4950 5000 5050 

Wavelength (A) Wavelength (A) 



Figure 1. Fits to the [OIII] emission lines profile of all the objects in our sample, with the exception of IRAS F13451+1232, IRAS 
F00188-0856 and IRAS F23327+2913. In the case of IRAS F13451+1232 the fit to the [OIII] profile can be found in Holt ct al. (2003) 
and therefore is not included in the figure. In the latter two cases the [OIII] emission lines have low equivalent widths and we used the 
Ha+[NII] complex to find one model that reproduces all emission lines. In each case the best fitting model is over-plotted (bold line) on 
the extracted spectrum (faint line). The different kinematic components are also plotted in the figure. The narrow, intermediate, broad 
and very broad (if present) components correspond to the the solid line, dashed line, dotted line and dot-dashed line. The exceptions 
arc IRAS F13305-1739, IRAS F15130-1958 and IRAS 23060+0505. In the first case the dashed line corresponds to one of the two 
broad components detected in this object, refereed as Bl in Table 2. In the case of IRAS F15130-1958, the solid line corresponds to the 
intermediate component referred as II in Table 2. Finally, in the case of IRAS 23060+0505, the dashed line corresponds to the narrow 
component refereed as N2 in Table 2. In addition, IRAS F14394+5332 shows a spectacular [OIII] profile that is not observed in all the 
other emission lines detected in its spectrum. See Notes on individual sources for details of the modelling for this and the other ULIRGs 
in our sample. 
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Figure 1. Continued 
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IRAS F23327+2913* IRAS F23389+0300 
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Wavelength (A) Wavelength (A) 

Figure 1. Continued 
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Table 3. 



Object H/3 NeVA3426 [OII]A3727 NcIIIA3869 [OIII]A5007 [OI]A6300 Ha A(Ha/H/3) [NII]A6583 [SII]AA6716,6731 

IRAS erg cm' 2 s _1 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) 



F23233+2817 e 


N 


8.10E-16 


0.45 


1.28 


0.42 


1.98 


0.33 


5.12 


0.69 


5.66 


2.37 




I 


8.86E-16 


0.88 


2.59 


1.26 


10.68 


0.63 


4.46 


1.12 


13.37 


1.41 




13 


1.58E-15 


1.10 


0.41 


1.14 


7.34 


0.40 


2.10 


0.33 


4.99 


1.04 


F23327+2913^ 


N 






0.14 




0.05 


0.08 


5.82E-15 


0.65 


0.47 






I 


3.02E-15 




0.91 




0.62 


0.26 


2.55 


0.14 


4.68 


0.75 


F23389+0300 6 


N/I 


4.12E-16 




3.33 


0.28 


2.28 


4.88 


7.87 


1.02 


4.71 


2.52 




VB 


9.88E-16 




1.91 


0.43 


3.63 


7.34 


8.35 


1.26 


25.88 


9.46 



Table 3. Table showing the fluxes of various emission lines detected in the optical spectra of the sources. Col (1): object name. Col (2): the label of the different components (Same as 
Col (2) in Table 2). Col (3): H/3 flux in erg s _1 cm -2 . From Col (4) to Col (9) and columns (11) and (12): ratios of the fluxes of the main emission lines to the H/3 line. The uncertainties 
in these ratio are typically < 15%, as estimated accounting for the 5% flux calibration uncertainty, and the uncertainties in the model fits themselves. However, there are a few particular 
cases, mainly the broadest component in some emission lines, or those kinematic components that make a small contribution to the flux, for which the uncertainty can be as high as 50% 
(for example, the 12 component of the [SII] AA6716,6731 emission lines in the case of IRASF15130-1958). Col (10): uncertainties in the H«/H/3 ratio (we discuss in the text a possible 
correlation between the broadness of the different kinematic components and their corresponding reddening. Therefore, it is particularly important to give the uncertainties associated 
with this ratio explicitly). 

a As a first approach during the modelling, [Oil] AA3726, 3729 was treated as a singlet rather than a doublet. However, there are cases for which accounting for the doublet nature of at 
least one of the different kinematic components is necessary to adequately model the [Oil] AA3726, 3729 emission lines. The fluxes presented in the table for these cases are are the sum 
of the two lines of the doublet. ^ 
b For these four galaxies it was not possible to find one model that adequately reproduces all the emission lines. In the cases of IRAS F13451+1232, F16156+0146 and F23389+0300 the 5r 
number of components is the same for all the emission lines, although with different widths and shifts. Therefore, it is still possible to use the H/3 line as a reference and give line fluxes _^ 
relative to that line. In the case of F14394+5332 the modelling results reveal the presence of a narrow component at rest frame found in all the emission lines, plus a number of different 
blucshiftcd components with different widths and shifts for the different emission lines. Therefore, to give an idea of the fluxes from the narrow component and the shifted components, g" 
we define two kinematic components in this case: narrow (N), corresponding to the narrow component found for all emission lines, and broad (B), which comprises the sum of all the 
kinematic components for each line that are not the narrow component. q 
c This is the only object with a strong broad-line Syl nucleus the in our sample. Due to the large contribution of the broad component from the broad line region (BLR) in the case of 
the Ha emission, it is not possible to constrain the flux in the broad, blueshifted component detected in [OIII] (referred as B in the table) for this particular emission line. Therefore, 5^ 
no flux value for the B component is presented in the table for this galaxy (sec appendix A for details) . 

d Two models adequately modelled the emission lines in the case of F17179+5444. One model comprises a narrow, an intermediate and a broad components, while the second model g 
includes two narrow components of the same width plus a broad one. The results shown in the table for this galaxy correspond to the first model. 

e The very broad (VB) kinematic component for these two ULIRGs is only detected in high ionization emission lines (i.e. not in the case of H/3) and therefore not included in the table. § 
f The H/3 emission line for this galaxy falls close to the region of the spectrum that is heavily affected by the dichroic. Therefore, the narrow component, which makes a relative small 
contribution for most of the emission lines, is not detected in the case of H/3. For this components line fluxes arc given relative to Ho. 2 
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3.2 The incidence of kinematically disturbed 
components in Seyfert and non-Sy-ULIRGs 

Kinematic disturbance of the emission line gas does not nec- 
essarily involve large line shifts, but can also be indicated by 
large line widths in components that are relatively unshifted. 
Therefore we define the term "kinematic disturbance" in 
this paper to encompass kinematic components that have 
intermediate or broad line widths (FWHM > 500 km s _1 ) 
and/or large blue shifts relative to the galaxy rest frame 
(AV < —150 km s" 1 ). By this definition, the overwhelm- 
ing majority (16/17 or 94%) of the Sy-ULIRGs for which 
we have been able to measure the emission line profiles in 
this study show evidence for kinematic disturbance in their 
nuclear regions. In contrast, only a small minority (3/16 or 
~19%) of the non-Sy-ULIRGs from Rodriguez Zaurin et al. 
(2009) present evidence for kinematic disturbance in the nu- 
clear 5 kpc aperturefl This difference is significant at the 
~3cr level, strongly suggesting that the kinematically dis- 
turbed components are driven by the AGN in most cases, 
rather than by the circum-nuclear starbursts. Significantly, 
the one Sy-ULIRG that shows no clear signs of nuclear kine- 
matic disturbance — F12112+0305 — is also the object that 
presents the weakest evidence for Seyfert-like AGN activ- 
ity based on its line ratios — falling close to th e AGN/HII 
compo site region in the diagnostic diagrams of lYuan et al.l 
|2010l ); it also has the lowest [OIII] emission line luminosity 
of all the Sy-ULIRGs in our sample (see Table 1). 

At first sight these result s may seem in contrast with 
those of ISoto fc Martini (|2012h and ISoto et all i|2012l ) who 
studied a sample of 39 ULIRGs and found evidence for dis- 
turbed kinematics in the majority of the objects in their 
sa mple, which is d o minat ed by non-Sy ULIRGs. However, 
the lSoto fc Martini (|2012l ) spectroscopic dataset has higher 
spectral and spatial resolution than ours. In addition, the 
sizes of the extraction apertures are different from those used 
in this paper. Therefore, it is likely that they detect some 
broad components th a t are not o bserved in our spe ctra. Fi- 
nally, ISoto fc Martini l|2012l) and ISoto et all (|2012h used a 
different modeling approach from the one we use here (i.e. 
a "simplest model approach"). As a result, their kinematic 
models for the nuclear emission lines in some of the sources 
in their sample comprise more than one kinematic compo- 
nent, regardless of whether the additional components are 
actually required to produce an adequate fit. 

3.3 The diversity of emission line profiles 

Figure 1 and Table 2 demonstrate a striking diversity 
in the emission line profiles of Sy-ULIRGs. The profiles 
range from the complex, multi-component profile visible 
in F14394+5332, through the case of IRAS F17044+6720, 
in which there is a single broad, blueshifted component 
along with a strong narrow component, to the case of 
F23389+0300, which has a strong, and extremely broad 
component that makes up a large proportion of the flux, yet 
is barely shifted relative to the narrow component. The min- 
imum number of Gaussian components required to fit the 

4 These three ULIRGs are F 14348-1447 and F23234+0946 and 
F15327+2340. 



individual emission lines ranges from 2 to 5, with 7 objects 
(~41%) requiring two, 7 objects (~41%) requiring three, 2 
objects (~12%) requiring four, and 1 object (F14394+5332) 
requiring five components for an adequate fit. However, it 
is important to note that, in the two cases - F23060+0505, 
F23233+2817 - where a fourth component is required to ad- 
equately model the emission lines, the relative flux contribu- 
tion of this fourth component is relatively small. Therefore, 
its FWHM and shift with respect to the narrow component 
are relatively unconstrained. 

In terms of determining the direction and magnitude 
of the velocity shift of any kinematically disturbed compo- 
nent, it is important to carefully determine the rest frame 
velocities of the ULIRG host galaxies (see Holt et al., 2003, 
2006, 2008). In the absence of other indicators (e.g. accu- 
rate stellar kinematics) the best way to do this is to com- 
pare the kinematics of the extended narrow components 
along the slit with the various kinematic components de- 
tected in the spatially-integrated emission line profile of the 
5 kpc aperture. For example, if the undisturbed gas shows 
a rotation curve across the nucleus, the component in the 
spatially-integrated emission line profile that corresponds to 
the galaxy rest frame will have a wavelength corresponding 
to the centroid of the blue and redshifted components on 
either side of the nucleus. In the majority of the cases in 
which we have spatially-resolved information on the narrow 
line kinematics, we find that the redmost narrow compo- 
nent in the spatially integrated profile coincides with the 
galaxy rest frame. However, in two cases — F17179+5444 
and F23060+0505 — a narrow and a relatively low-FWHM 
intermediate component (FWHM = 515 ± 33 km s _1 ), or 
two narrow components, are present in the nuclear line pro- 
files, each corresponding to the velocity shift of the spatially 
resolved gas on one or other side of the nucleus. In these 
cases, the observed narrow line splitting is likely to repre- 
sent unresolved rotation, or to trace a large-scale b ipolar 
outflow in the gas (see, for example iHolt et ail . 120081 ). In 
these two cases we used the average of the wavelengths of 
the two narrow components to define the galaxy rest frame. 

With the exceptions of F17179+5444 and 
F23060+0505, the velocity shifts listed in Table 2 are 
defined relative to the redmost narrow component fitted in 
the 5kpc aperture spectrum. It is clear that the broadest 
kinematic components in each object are significantly 
blueshifted relative to the galaxy rest frames in all cases 
except F13428+5608 and F23389+0300. In the case that 
the emission line clouds suffer extinction by a general 
field of dust in the nuclear regions of the galaxies, rather 
than by dust in the clouds themselves, this implies that in 
most cases the gas is undergoing a radial outflow from the 
nucleus. 

Finally, we can define the flux weighted mean line width 
and velocity shift of the disturbed kinematic components as: 

FWHM, = ^,FWH Mi) 
where Fi, FWHM; and AVi are the fluxes, FWHM and 
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Figure 2. Flux weighted mean FWHM (FWHM W ) and shifts 
(AV W ) plotted against the [OIII]A5007 luminosities (L[oiiij) of 
the sources. 

velocity shifts of all the intermediate and broad kinematic 
components with respect to the rest frame. The results are 
shown in Table [4] For our sample of 16 ULIRGs, we find 
values in the range of 797 < FWHM W < 2223 km s" 1 with 
a mean and median values of 1253 km s _1 and 1336 km 
s" 1 respectively. In the case of the AVW we find vaues in 
the ranee -942 < AU W < 47 km s _1 with mean and median 
values of -406 km s^ 1 and -333 km s _1 respectively. 

Figure [2] shows FWHM W and AVW plotted against the 
total [OIII]A5007 emission line luminosity (I/[oni]) from col- 
umn 6 of Table 1. This figure is designed to investigate any 
possible correlation between the properties of the outflows 
and the AGN power, as indicated by the I/[oni] luminosity, 
which is often co nsidered to be a reliable AGN bo l omet- 
ric in dicator (e.g. D icken et all |2010| ; lLaMassa et all 120091 . 
2010). However, we do not find any clear correlation be- 
tween Lroiii] and either the shifts or the widths of the out- 
flows observed at optical wavelengths, for the galaxies in 
our sample. Our results suggest that the kinematic prop- 
erties of the AGN-induced outflows present in Sy-ULIRGs 
are not strongly dependent in the instantaneous AGN radia- 
tive power as indicated by the [OIII] luminosity. A possible 



Figure 3. Same as Figure[2]but with the 1.4 GHz monochromatic 
radio luminosities (Li.4GHz)- The dashed line indicates the maxi- 
mum 1.4 GHz luminosity for the galaxies in our sample calculated 
using their Lir along with the Condon et al. (1991) IR-radio lumi- 
nosity relation for starburst galaxies. Four objects show an excess 
of radio emission for their corresponding Lir, which reveals the 
presence of a strong AGN component. These objects are indicated 
in the Figure using the first 5 digits of their designation in the 
IRAS Faint Source Catalogue Database. 



reason for this lack of a correlation is that the [OIII] lumi- 
nosity represents the AGN power on the timescale of the 
light crossing time of the narrow line region (iNLR < 10 4 yr), 
but the AGN activity may be highly variable within the 
~ 10 8 yr timescale of the final stages of gas-rich mergers 
represented by the ULIRGs. In this case, the instantaneous 
[OIII] luminosity might not provide a good indication of 
average AGN power over the entire activity cycle. Alterna- 
tively, because of the substantial extinction on a kpc-scale 
in the NLRs of ULIRGs, the [OIII] luminosity may not be 
as reliable an indicator of the AGN bolometric luminosity as 
it is in other types of active galaxies. Indeed, many ULIRGs 
show co mpelling evidence for AGN activity at mid-IR wave- 
lengths jVeilleux et all [2009). but no sign of such activity 
at optical wavelengths. 

Unfortunately, estimating the reddening is challenging 
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Table 4. Weighted mean FWHM and velocity shifts (FWHM V 
and AV W ) for the objects in our sample. The final column gives 
the estimated radii of the [OIII] outflow regions. For most of the 
objects these have been estimated by fitting a Gaussian to the core 
of the [OIII]AA5007,4959 spatial profile along the slit, with the 
radius estimated as -R[om] = FW H M(kpc)/2. However, in the 
case of F13451+1232 the radius was estimated fr om the narrow- 
band images published inlBatchcldo r et al. | J2007| V 



based on our existing data. Ha is in a blend, and there are 
potentially degeneracies involved in the fits that may af- 
fect the accuracy of the line ratio measurements, especially 
in cases with highly complex, multiple-component line pro- 
files. Nonetheless, we used the Ha/H/3 ratios in Table 3 to 
calculate a "crude" estimate of the de-redenned [OIII] lu- 
minosities. We then used these de-reddened luminosities to 
reproduce Figure 2 and obtained identical results, i.e. no 
trends and correlations are found. Furthermore, we carried 
out t he same excerc i se but using the [OIV]25.89 luminosities 
from IVeilleux et all (2009), which are a more reliable AGN 
luminosity indicator and are available for 9 ULIRGs in our 
sample. The results obtained are identical to those of the 
previous attempt, i.e., no trends and correlations are found. 



3.4 The spatial scales of the kinematically 
disturbed regions 

Our analysis of the emission line profiles has concentrated 
on spatially integrated, 5 kpc diameter aperture, data. 
However, in order to properly quantify the outflows it 
is important to attempt to estimate their spatial extents 
within the apertures. Spatially resolved, kpc-scale emission 
line outflows are found in many Seyfert galaxies, and di- 
rectly detected in one low-z S y-ULIRG (Mrk23 1: see R upke 
fc Veilleux 20 1 1). M oreover, lAlexander et al.l (|2010l l and 
lHarrison et al.l (|2012j 'l have recently found evidence for out- 
flows on a scale of 4 - 15 kpc in a sample of high-z ULIRGs. 
In contrast, HST emission line imaging of the outflows in 
two radio-loud ULIRGs - PKS1345+12 and PKS1549-79 - 
shows that they are significantly more compact, with ra - 
dial extents of only ~100 - 200 pc l|Batcheldor et alllioOTt) . 



Therefore, the radial extents of the AGN-induced outflows 
in the Sy-ULIRGs in the local Universe are highly uncertain. 

Based on visual inspection of our 2D long-slit data, we 
find that the dominant broad/intermediate kinematic com- 
ponents are generally confined to the central few pixels - on 
a similar scale to the seeing FWHM estimated by the DIMM 
seeing monitor at the time of the observations. More quanti- 
tatively, we have extracted the continuum-subtracted spatial 
profiles of the [OIII] emission from the 2D frames and fitted 
Gaussians to the central cores of the profiles. We show the 
results in column 4 of Table 4, quantified as FWHM/2 in 
units of kpc. We find that the [OIII] FWHM estimates are 
comparable with the DIMM seeing estimates at the times 
of the observations for ~56% of the objects in our sample. 
Moreover, in some objects the spatial profiles may be dom- 
inated by spatially resolved narrow components. Therefore, 
in most cases the numbers presented in column 4 of Table 
4 are likely to represent upper limits on the true spatial 
extents of the near-nuclear warm outflows. 

Overall, based on our measurements of the spatial pro- 
files of the [OIII] emission in our sample of Sy-ULIRGs, and 
the direct HST imaging estimates of the sc ales of the outflow 
regio ns in PKS1345+12 and PKS1549-79 (|Batcheldor et all 
120071 ) . it is likely that the radial extents of the near- nuclear 
outflows fall in the range 0.1 < r < 3.5 kpc. 



4 DISCUSSION 

4.1 Are the outflows jet-driven? 

ISpoon fc Holt] ( 2009) found evidence for a correlation be- 
tween the width of the mid-IR [Nelll] line at 15.56 /im and 
the 1.4 GHz radio luminosity (I/i.4GHz), based on Spitzer 
IRS observations of a combined sample of ULIRGs and lo- 
cal Seyfert galaxies. The correlation is in the sense that some 
of the objects with the highest radio powers have the most 
extreme emission line kinematics, suggesting that the out- 
flows might be jet-driven. Further evidence for jet-driven 
outflows in ULIRGs is provided by HST images which show 
that the outflow regions have similar scales to the relativis- 
tic radio jets in two of the m ost extreme radio-loud U LIRGs 
with high velocity outflows (Batchcld or et ail . l2007t ). 

In Figure [3] we show the weighted mean velocity shifts 
and FWHM plotted against the radio powers listed in col- 
umn 7 of Table 1. Note that an important caveat with in- 
terpreting this plot in terms of jet-driven outflows is that 
the radio emission does not necessarily have an AGN ori- 
gin in the ULIRGs, but may instead be associated with the 
prodigious starburst activity. The vertical dashed line in [3] 
show the maximum radio power expected for the starburst- 
related radio emission in the Sy-ULIRGs, as predicted on 
from their measured infrared luminosities and the well- 
known correlation between radio power and inf rared lumi- 
nosity for starburst galaxies |Condon et al.l . IT991). From this 
it is clear that only four of the Sy-ULIRGs - F13305-1739, 
F13451+1232, F17179+5444, F23389+03000 - can be con- 
sidered to have genuinely radio-loud AGN with L\ 4GHz > 
10 24 ' 3 W Hz" 1 . 

On the basis of Figure there is no evidence for a 
strong correlation between radio power and emission line 
kinematics. Although it is certainly true that the two most 
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radio-powerful objects in our sample - F13451+1232 and 
F23389+03000 - also display some of the most extreme emis- 
sion line kinematics, the other two radio-loud Sy-ULIRGs 
show relatively quiescent kinematics. Moreover, the radio- 
quiet majority of objects in our sample (Z/i.4GHz < 10 24 W 
Hz - ) encompass the same range of kinematic disturbance 
as the radio-loud objects, and some radio-quiet objects 
show emission line kinematics as extreme as those found 
in the radio-loud objects, with [OIII] emission line profiles 
dominated by broad, shifted components (e.g. F01004-2237, 
F14394+5332, F15130-1958, F15462-0450). Therefore we ar- 
gue that, while the emission line outflows can be driven 
by the jets in some Sy-ULIRGs, especially in cases where 
powerful relativistic jets have the same spati al scales as 
the emission line struc t ures (e.g. F13451+1232 iHolt et all . 
120031 ; iBatcheldor et all 120071 ), in most Sy-ULIRGs there is 
no clear evidence for such jet-driven outflows. 



4.2 Ionization mechanisms 

4-2.1 AGN and starburst photoionization 

It is well known that starburst and AGN activity can 
co-exist in the nuclear regions of ULIRGs, both energy 
sources potentially contributing to th e ionization the sur- 

1 19981: iFarrah et all 120071 : 



rounding gas (e.g. iGenzel et al 



lArmus et all . 120071 ; IVeilleux et al 



2009). Therefore it is in- 



teresting to determine whether there there is correlation be- 
tween the emission line kinematics and the nature of the 
ionizing source(s). 

Figure [4] shows the standard, optica l BPT diagnostic 
line ratio diagrams (Bal dwin et all Il98ll ). The lines drawn 
in the di agrams correspond t o the the optical classification 
scheme of lKewlev et alj j2006) . Overplotted on the figure are 
the results from our kinematic study. Blue squares, green 
triangles and red circles correspond to the narrow, inter- 
mediate and broad components respectively. The error bars 
have been calculated by combining the 5% uncertainty of the 
relative flux calibration in quadrature with the uncertainty 
associated with the model fits. 

It is notable the line ratios obtained for most, if not all, 
the ULIRGs in our sample are consistent with a Sy2 clas- 
sification. However, the figure demonstrates a clear trend 
with emission line kinematics: whereas the line ratios for the 
broad components (red circles) in most cases fall squarely in 
Seyfert part of the diagram, suggesting AGN photoioniza- 
tion as the dominant ionization mechanism, the line ratios 
for the narrow components (blue squares) measured in many 
objects fall in, or close to, the transition zone between HII 
region and Seyfert classifications - more consistent with a 
combination of the AGN and stellar photoionization; the 
ratios measured for the intermediate components (green tri- 
angles) tend to fall in between those of the narrow and the 
broad. 

We emphasise that the trends in the emission line ra- 
tios with linewidth do not necessarily provide evidence for an 
ionization stratification in the AGN outflow regions them- 
selves. Rather, they are naturally explained in terms of the 
varying contribution of the starburst activity to the ioniza- 
tion of the gas associated with each of the different kinematic 
components. For example, it is likely that starburst activity 
plays an important role in ionizing the gas associated with 



the narrow components in many of the objects. The impor- 
tance of the starburst activity as an ionizing source grad- 
ually decreases for the intermediate component. Finally, in 
the case of the broadest kinematic component AGN activity 
is the dominant source of ionization. Note that the associ- 
ation between the ionization state of the gas and the line 
width provides further evidence that the broader (generally 
blueshifted) components are driven by the AGN rather than 
by the starburst components. 

With the aim of better understanding the nature of ion- 
ization mechanisms resp onsible for the measu red line ratios, 
Figure [4] also shows the I Groves et alj (l2004aT ) dusty, radia- 
tion pressure-dominated photoionization models for NLG in 
AGN. For these models we assume a hydrogen density n# 
= 1000 cm -3 (the only hydrogen density used in the Groves 
et al. 2004a paper) and t wice solar abundance (2Zq). This 
later value was found by I Groves et al.l (I2004br ) to provide 
a better mat c h to th e line ratios measured in Seyfert NLR 
iGroves et al.l (|2004bh . A grid of models for various values of 
the ionizing parameter (Uo — 0, -1, -2, -3, -4) and ionizing 
continuum SED power-law indicies (F„ <x v a , a = -1.2, -1.4, 
-1.7, -2.0) is shown in the figure. 

At first sight it is apparent that, with few exceptions, 
the line ratios of the intermediate and the broad kinematic 
components fall within the grids of models in the case of 
the diagram involving [OI]/Ha. However, these models fail 
to reproduce most of the line ratios for such components 
in the case of the [Nil] /Ha diagram (and to a lesser ex- 
tent in the case of the [SII]/Hq diagram). A possible ex- 
planation for this finding is that the emission line regions 
are more metal rich than we have assumed. In this context 
Figure[S]shows the [OIII]5007/H/3 vs [NII]6583/Ha diagram 
assuming njj = 1000 cm -3 and 4Zq (the other two diagrams 
remain relatively unchanged after increasing the metallicity 
and therefore are not shown in the figure). It is clear from the 
figure that a 4Zq abundance certainly improves the match 
between the models and our data. Therefore, the line ratios 
derived for the Sy-ULIRGs in our sample are consistent with 
gas of super-solar abundances that has been photoionized 
by a combination of starburst and AGN components, with 
the AGN photoionization dominating for the broad/shifted 
kinematic components. 



4-2.2 Shock ionization 

As well as photoionization by AGN and starburst activ- 
ity, it is possible that line ratios plotted in Figure [4] can 
also be ex plained in terms of i onization by fast radia- 
tive shocks dDopita fc Sutherland! Il995l . 1 19961 : IGroves et all . 
l2004al : I Allen et all 120081 '). Theshock models predict a se- 
ries of line ratios for a range of magnetic field (B), electron 
densities (n e ), abundances and shock velocities (v s ). In ad- 
dition, the predicted line ratios depend on the geometry of 
the shock, i.e. the presence of a photoionized precursor of 
the shock f see I Allen et al.. 2008, for a de tail discussion). 

Figure[5]shows the Allen et al.l (120081 ) models with solar 
metallicity (Zq) and a pre-shock density of 100 cm" 3 . Based 
on the results in the previous section, it might be more ap- 
propriate to use models with hig her metallici t ies. H owever, 
the highest metallicity used in the I Allen et all l|2008h is 2Z Q , 
and for this value only n e values of 1 and 1000 are consid- 
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Figure 4. Standard optical diagnostic diagrams showing the c l assific ation scheme by iKewlev et al. The solid curves are the 

theoretical "maximum starburst line" derived by IKewlev et al. I 1120011 ) as an upper limit for star-form ing galaxies and the em pirical 
boundary lines between Seyfert 2 galaxies and LINERs. Th e dashed curve on th e [Nil] diagram is the iKauffmann et al.l J2003f) scmi- 
empirical lower boundary for the star forming galaxies. The I Groves et al . (2001a) grids of dusty, radiation pressure-dominated models 
are also plotted in the Figure. These grids have been generated assuming twice solar metallicity (2Zq) and hydrogen density of njj = 
1000 cm -3 . Gridlines corresponding to five values of ionizing parameter (Uo = 0,-1,-2,-3,-4) and four values of power law index (F„ oc 
v a , a = -1. 2, -1.4, -1.7, -2.0) are shown in the figure. To help the reader follow the gridlines these are grey-coded from "light-grey" to 
black, with light-grey and black corresponding to the lowest and highes values of Uq and a respectively. Over-plotted on the diagrams are 
the results of our kinematic study. Blue squares are the line ratios corresponding to the narrow components, green triangles correspond 
to the intermediate component and red circles represent the broadest kinematic component. 



ered. Therefore, we decided to use Zq, models since they 
sample the parameter space better. 

One of the largest uncertainties in these models is pre- 
shock density. Assuming that we measure the density of the 
compressed post-shock gas in the kinematically disturbed 
emission line components that we detect in the ULIRGs in 
our sample, and the measured electron densities of these 



components are 1000 cm -3 (see section 4.3), then the pre- 
shock densities could be as low as 10 - 100 cm -3 , since 
the compression fa ctor in the cooled, pos t -shoc k gas can be 
high (~10-100, see lDopita fc Sutherland! Il995h . Therefore, 
we have assumed a p r e-shock density of 1 00 cm " 3 . Finally, 
iDopita fc Sutherland! l|l995l ); lAllen et al.1 (120081 ) found that 
the shock+precursor models produce a better fit to the line 
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Dopi t a and collaborators ( e.g. Dopita fc Sutherland! . Il995l . 
1 19961 ; iGroves et all l2004bl ; I Allen et all . 120081 ), it is not al- 
ways possible to decisively distinguish between shock ioniza- 
tion and AGN photoionization using only the stronger opti- 
cal emission lines included in the diagnostic plots of Figure 



-1.0 -0.5 0.0 
log([NII]6583/Ha) 

Figure 5. Same as Figure [4] but only for the [NII]6583/Hq di- 
agram and for 4Zq. The other two diagrams remain relatively 
unchanged by the increase in metallicity and therefore are not 
shown in the figure. 



ratios measured for their samples of Sy2 and luminous IR 
galaxies. For this reason, we decided to use such models for 
the work presented in this section. The potential effects of 
changing these assumptions is discussed later in the section. 
The gridlines correspond to 3 values of B from 1-100 [iG (1, 
10 and 100 mG), which a re appropriate for starburst galaxies 
l|Thompson et all bOOrj ). and 9 values of v s (200, 300, 400, 
500, 600, 700, 800, 900 and 1000 km s" 1 ). The direction in 
which B increases is indicated in the figure with an arrow, 
while v s increases from bottom to top with 3 values (200, 
400 and 1000 km s" 1 ) indicated in the figure. 

As seen in the figure, most of the measured line ratio 
values do not fall within the region covered by the grids 
of shock models; the lack of agreement between the models 
and the data is particularly apparent in the diagnostic dia- 
gram that involves the [OI]/H« ratio. However, the models 
have a relatively large number of free parameters (B, v 3 , n e 
and abundance) plus the additional possibility of including 
or otherwise the shock precursor component. Therefore, it is 
possible to cover wide ranges of line ratio values using differ- 
ent assumptions. For example, using twice solar abundance 
models with a pre-shock density of 1000 cm" 3 instead of 100 
cm " 3 substantially ex tends the gridlines in the 3 diagrams 
( see I Allen et~aill2008l . for details), providing a better match 
to the measured line ratios. We note that the line ratio val- 
ues measured f or some of the bro ad kinematic components 
fall outside the I Allen et al.l (|2008T ) grids of shock models, or 
are only reproduced by the most extreme set of parameter 
values (e.g. n e = 1000 cm" 3 and B = 1000 ^G). 

Overall, we find that the line ratios measured for our 
sample of Sy-ULIRGs are, in general, better reproduced by 
photoionization models that include a combination AGN 
and starburst photoionization and in which the emitting 
gas highly enriched in metals. However, as pointed out by 



4.3 Physical conditions and reddening 

Potentially, further clues to the nature of the emission line 
outflows may be obtained by considering the reddening 
and densities o f the emissi on line components. For exam- 
ple, (|Holt et al.l . l2003l . l20lih found evidence that the broad- 
est and most highly shifted emission line components in 
F13451+1232 (PKS1345+12) are also the most highly red- 
dened, suggesting that the outflows are stratified and decel- 
erating, with the fastest components associated with most 
highly reddened regions closest to the driving AGN. Sim- 
ilarly, if the most highly disturbed kinematic components 
are truly associated with the regions closest to the AGN, 
we might expect them to show evidence for higher densities 
than the less kinematically disturbed components. 

On the outset, it is important to emphasize that esti- 
mating the reddening and densities for the broad, shifted 
kinematic components is challenging based on our exist- 
ing data, because some of the key emission lines (e.g. Ha, 
[SII] AA6717, 6731) are in blends, and there are potentially 
degeneracies involved in the fits that may affect the accu- 
racy of the line ratio measurements, especially in cases with 
highly complex, multiple-component line profiles (see dis- 
cussion in section 3.1). 

To search for reddening gradients, we use the Ha/H/3 
ratio, which, for typical NLR conditions and AGN photoion- 
izing continuum shapes, is expected to have a value in the 
range 2.8 < Ha/Hj3 < 3.1 in the absenc e of reddening 
( Gaskell fc Ferlandl . 11984k lOsterbrockl . Il989h . Column 9 in 
Table 3 gives the Ha/Hfi line ratio values for the vari- 
ous kinematic components of the galaxies in our sample. 
At first sight, the modelling results appear to show evi- 
dence for a trend of increasing reddening with linewidth in 
four of the ULIRG Seyferts (F13451+1332, F14349+5332, 
F17179+544, F23060+0505). However, two objects appar- 
ently show the reverse trend (F00188-0856, F15462-0450), 
and in most cases there is no clear trend. Clearly, given the 
large uncertainties in the measured ratios — to some ex- 
tent reflecting the degeneracies in the fits — the results are 
inconclusive. 

The results are no more decisive when it comes to 
estimating the densities using the [SII] (6717/6731) ratio 
(see the discussion of individual objects in the Appendix). 
While it has proved possible to derive accurate densi- 
ties for the narrow emission line components in some ob- 
jects (typical electron densities: 2.4 x 10 2 < n e < 3.6 x 
10 3 cm" 3 ), we have failed to measure accurate densities for 
the broader emission line components in any object. How- 
ever, in the cases of F13305-1739, F13428+5608, F15462- 
0405 and F23233+2817, it has been possible to derive lower 
limits on the electron densities of the broad components of 
9xl0 3 , 6xl0 3 , 4xl0 3 and 4xl0 3 cm" 3 respectively. 
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n e =100 crrT 3 ,Z n e =100 cm~ 3 ,Z 




log([NII]6583/Ha) log([0I]6300/Ha) 

n e =100 crrT 3 ,Z 



1 .5 




log([SII]6716,6731/Ha) 

Figure 6. Same as Figure [4] but showing the Al len et al. grids of shock-ionization models. These grids have been generated 

assuming a shock+precursor model with solar mctallicity and prcshock density of n e = 10 cm -3 . Gridlines corresponding to 3 values of 
magnetic field (B = 1, 10 and 100 fiG) and 9 shock velocity values (v s = 200, 300, 400, 500, 600, 700, 800, 900 and 1000 km s^ 1 ) are 
shown in the figure. The increasing direction of B is indicated with an arrow. In the case of the [SII]/Ha diagram B increases first to the 
right and then to left which is indicated with a curved arrow. v s increases always from bottom to top with 3 values (200, 400 and 1000 
km s — 1 ). As in the case of Figure[4] the gridlines these are grey-coded from "light-grey" to black, with light-grey and black corresponding 
to the lowest and highest values of B and v s respectively. 



4.4 Mass outflow rates and energetics 



AGN feedback is now routinely incorporated into numerical 
simulations of the hierarchical growth of galaxies through 
major galaxy mergers. To reproduce the correlations be- 
tween the masses of the black holes and the bulge prop- 
erties (e.g. ISilk fc Reesl . Il998l : Idi Matteo etafl . 120051 ) some 
of these models require that a relatively large fraction of 
the available accretion power of the quasars is thermally 



coup l ed to the circumnuclea r ISM (~ 5 - 10%, e.g. iFabianl . 
Il999l ; ldi Matteo et al.l . l2005l . but see Hopkins et al. 2010). 

As described in the introduction, ULIRGs represent just 
the situation modelled in many of the merger simulations 
that include AGN feedback. In order to investigate whether 
our observational results are consistent with those of the 
merger simulations, it is important to quantify the mass 
outflow rates and the kinetic powers of the outflows. 

For a spherical outflow, the mass of the ionized gas in 
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the outflow is given by: 



M = n e m p Ve 



(3) 

where n e is the electron density, m p the mass of the pro- 
ton, V is the total vo l ume and e is the gas filling factor 
( Oste rbrock fc Ferlandl 12006). Therefore, for a steady state 
outflow, the mass outflow rate (M) is: 



M = n e m p VoutAe 



(4) 



where v ut is the velocity of the outflow and A is surface area 
of the outflowing region. Also, filling factor e is related to 
the H/3 luminosity, the electron density n e and the emitting 
volume V by: 



L(H/3) 



(•») 



hence, the mass outflow rate is related to the H/? luminosity 
by: 



M 



L(R/3)m p v utA 



v J J 
*H/3 



For a spherical outflow geometry this reduces to: 



M = 



3L(Ul3)m p v out 



(6) 



(7) 



where r is the radius of the spherical volume. In addition, 
the kinetic power of the outflow (E) is related to the veloc- 
ity dispersion (a w FWHM/2.355), mass outflow rate and 
outflow velocity by: 



E=y(V& t + 3a 2 ) 



(8) 



The main uncertainty in calculating the mass outflow 
rates, and hence the kinetic powers of the outflows, is related 
the electron density. Unfortunately, as we saw in the last 
section, due to degeneracies in the fits to the [SII] lines, the 
densities are relatively unconstrained for the broad, shifted 
components in most of the objects in our sample. In this 
case, the densities could range from the relatively low val- 
ues typically determined for the spatially-resolved NLR in 
some objects (n c ~ 100 cm" 3 : Taylor et al. 2003, Robin- 
son et al. 2000) up to the much higher densities directly 
measured for the outflow components in F13451+1232 us- 
ing th e transauroral lin e ratios (2 x 10 4 < n c < 5 x 10 5 
cm - 3 : lHolt et aill201ll ) - a range of almost four orders of 
magnitude. Other sources of uncertainty include the H/3 lu- 
minosities (Lh/s), which are affected by uncertain reddening 
corrections (see the last section), and the uncertain radii (r) 
of the outflows (see section 3.4). Each of these adds a least 
a further order of magnitude uncertainty when calculating 
the properties of the outflows. 

In order to gain an impression of the range of possible 
mass outflow rates and kinetic powers for the warm outflows 
in the Sy-ULIRGs, we have calculated upper and lower limits 
on these quantities, based on the emission line kinematics 
from Table 2 and assuming that AV = V ou t for each of the 
kinematic components. To derive the lower l i mits w e have 
used the high density estimated bv lHolt et all (|201lT ) for the 
broadest component in F13451+1232 (n e = 5 x 10 5 cm" 3 ), 
the upper radius limits given in column 4 in Table [4] and 
H/3 luminosities that have not been corrected for reddening; 
whereas for the upper limits we have assumed a density of 



100 cm" 3 , a radius of 0.1 kpc, and the reddening-corrected 
H/3 luminosities. Note that low densities, small radii and 
high H/3 luminosities favour highly massive and energetic 
outflows. In the case of F13451+1232, the density, reddening 
corrected H/3 luminosity, and radius of the outflow are know 
relativ ely accurately f r om th e studies of iBatcheldor et all 
(|2007h and iHolt etafl (|201lh : for this object we have used 
the results from the latter studies to calculate the outflow 
properties. 

As expected, with the exception of F13451+1232, the 
estimates of the outflow properties are highly uncertain, 
and encompass several orders of magnitude in both mass 
outflow rate and kinematic power for individual kinematic 
components in individual objects. Across our sample we 
find upper limiting mass outflow rates in the range 1 < 
M < 1.5 x 10 3 Mq yr _1 and upper limiting kinetic pow- 
ers in the range 5 x 10 41 < E < 2 x 10 45 erg s" 1 . To 
put these numbers into context, we can compare them with 
the bolometric radiative powers of the AGN, which provide 
an indication of the available accretion power of the ma- 
terial being accreted onto the super-massive black holes. 
We calculate the bolometric luminosities of the individual 
AGN by multiplying the infrared luminosities from Table 1 
by the fractional contributions of the AGN to the infrared 
light f rom mid-IR sp e ctrosc opic study of the QUEST sam- 
ple bv lVeilleux et all l|2009l ff obtaining values in the range 
2 x 10 45 < Lb§l < 6 x 10 45 erg s" 1 . Therefore, considering 
the extreme high end of the upper limits, the warm outflows 
in Sy-ULIRGs can have kinetic powers that are compara- 
ble with the radiative luminosities of the AGN, However, 
at the lower end of the range of upper limits - includ- 
ing F1345 1 +123 2 with its more precise E estimates (see 
IHolt et al l . l201lt) - the kinetic powers of the outflows are 
only a small fraction of the bolometric luminosities (<0.5%) 
— well below the l evels required by some of the models (e.g. 
ISilk fc Reed . 1 19981 : iFabianl . Il999l : Idi Matteo et~aLl , l2005h , but 
perhaps cons istent with the two stag e feedback mechanism 
described bv lHopkins fc Elvisl (|201CJ ). 

The discrepancy with the models grows even stronger 
if we consider the range of lower limiting mass outflow rates 
and kinetic powers encompassed by our sample: 4 x 10~ 5 < 
M < 4 x 10" 3 M Q yr" 1 and 4 x 10 35 < E < 2 x 10 39 erg 
s _1 . If these values proved correct, the AGN outflows would 
have a negligible impact on the overall evolution of their 
host galaxies. 



4.5 Comparison with the emission line kinematics 
of other samples of Seyfert galaxies and other 
types of AGNs 

The hydrodynamical models predict that the AGN-induced 
outflows should be particularly strong in the final stages of 
the type of major, gas rich mergers represented by ULIRGs. 
Therefore it is interesting to consider whether the degree of 
kinematic disturbance in ULIRGs with optical Seyfert nuclei 
is unusual compared with samples of non-ULIRG AGN. 



5 Not e that, for objects in our sample that were not included 
in the lVeilleux et al] (2009) QUEST sample, we assume that the 
AGN contribute 50% of the infrared light - a typical value for the 
QUEST sample objects. 
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Sample Selection criteria Redshift range L[oiii] range (W) Resolution (km s 1 ) 



Sy-ULIRGs (this paper) 


S 6 o >1 Jy, h IR > 10 12 L Q , 
emission line ratios 


z < 0.175 


2 x 10 33 - 


- 4 x 10 35 


~230 


Heckman et al. (19811 


representative of 
local Seyferts 


z<0.076 


3.5xl0 32 


- 1.7xl0 35 


~130 


Whittle (1985) 


representative of 
local Seyferts 


z< 0.066 


2.3xl0 32 


- 3.3xl0 35 


40-80 


Veilleux (1991a) 


F[OIII] > 1x10" 13 
erg s _1 cm" 2 


z<0.028 


5.8xl0 32 


- 1.8xl0 35 


~10 


Nelson & Whittle (1995) 


representative of 


z< 0.043 


2.1xl0 31 


- 1.7xl0 35 


80-230 



local Seyferts 

Table 5. Details of the properties of th e local Seyfert compari s on samples. The [OIII1A500 7 flux values for the sources within these 
comparison samples ha v e been taken from:|Pah ari fc De Robertid dl988 h.lGu et al.l ll2006h.lMo iistakas fc Kennicuttl j2006l) , iMassaro et al.l 
ll2009l ) lTremblav et all J2009h . iGonzalez-Marti'n et al.l ||2009|) , lLaMassa et al.l d2010h . iGreene et al.l j201Cf) and Marvin Rose (private 
communication) . 



4-5.1 Comparison with previous kinematic studies of the 
NLR in Seyfert galaxies 

For our first c ompar iso n, we use the [ Heckman et "all 
Jl98ll)IWhittlel j 19851 ). iNelson fc Whittle! j 19951 ) and 
IVeilleuxl (|l991aT ) studies of the narrow emission ine kine- 
matics in samples of local Seyfert and radio galaxies. 
In general, although these samples are likely to be rep- 
resentative of nearby Seyfert galaxies, their selection 
criteria are not well defined. The exception is the IVeilleuxl 
l|l991af ) sample, for which the objects have S> -25° and 
F([OII Il <5007) > lxl0 _13 erg s" 1 cm" 2 (see IVeilleuxl . 
Il991al . for a discussion of the potential biases due to their 
selection criteria). The general properties of these samples 
are compared with those of the Sy-ULIRGs in Table [5] 
We note that, although these samples encompass a similar 
range of [OIII] emission line luminosity to our Sy-ULIRG 
sample, they tend to contain a higher proportion of lower 
luminosity objects. 

The previous studies of nearby Seyfert galaxies investi- 
gated the kinematics of the NLR using a parameter scheme 
based on cuts across the profiles at different heights. We will 
concentrate here in the so called asy mmetry (^4720) and kur - 
tosi s (Kl) parame ters as defined by I Heckman et all (|l98ll ) 
and lWhittld (|l985l ) as follows: 

, rnn WL2Q - WR20 

AI2 ° = WL20 + WR20 (9) 
where WL2Q and WR20 are the widths to the right and left 
of the line center defined at the 80% intensity level (C80), 
and: 

W50 



Kl = 1.524 x 



W20 



(10) 



where W50 and W20 are the widths at 50% and 20% inten- 
sity level (i.e W50 = FWHM of the line). The kurtosis Kl 
parameter is defined such that for a Gaussian Kl = 1.0. 
Line profiles with Kl values lower or higher than 1 are 
called platykurtic ("peaky") or leptokurtic ("stubby") re- 
spectively. 

Table E] shows the values of AT20, Kl, W20 and W50 
for the Sy-ULIRGs in our sample. Since the previous studies 
considered here concentrate mostly on the [OIII]A5007 emis- 
sion line, the values in the table have been estimated using 
that line. When more than one kinematic narrow compo- 
nent is clearly visible in the emission line profile of a source 



Object 


W20 


W50 


AI20 


Kl 


IRAS 


(km s" 1 ) 


(km s" 1 ) 






(1) 


(2) 


(3) 


(4) 


(5) 


F00188-0856 


1207±32 


509±56 


0.24±0.05 


0.64±0.10 


F01004-2237 


1779±30 


249±42 


0.71±0.03 


0.21±0.09 


F12072-0444 


1002±30 


544±36 


0.31±0.04 


0.83±0.14 


F12112+0305 


743±18 


388±36 


-0.20±0.06 


0.79±0.09 


F13305-1739 


1786±31 


1072±29 


0.11±0.02 


0.91±0.05 


F13428+5608 


949±36 


540±36 


0.00±0.04 


0.88±0.05 


F13451+1232 


706±31 


1553±36 


0.34±0.02 


0.68±0.03 


F14394+5332 


2091±42 


1488±30 


0.65±0.04 


1.08±0.05 


F15130-1958 


2029±30 


1156±42 


0.28±0.05 


0.87±0.14 


F15462-0450 


2208±46 


1659±40 


0.68±0.05 


1.14±0.09 


F16156+0146 


989±29 


297±30 


0.36±0.03 


0.46±0.13 


F17044+6720 


640±21 


296±29 


0.04±0.03 


0.70±0.03 


F17179+5444 


1202±29 


645±29 


0.07±0.03 


0.82±0.06 


F23060+0505 


1315±32 


721±30 


0.38±0.03 


0.83±0.10 


F23233+2817 


1459±34 


834±29 


0.31±0.05 


0.87±0.14 


F23327+2913 


1151±40 


639±40 


0.27±0.02 


0.84±0.11 


F23389+0300 


2294±36 


571±29 


0.07±0.01 


0.38±0.05 



Table 6. The emission line parameters for the Sy-ULIRSGs in our 
sample. Col (1): galaxy name. Col (2): width at 20% of the peak 
intensity level. Col (5): width at 50% of the peak intensity level. 
Col (4): the asymmetry yl/20. Col (5): the kurtosis parameter Kl. 
All these parameter values have been quadratically corrected for 
the instrumental profile. 



(e.g. IRAS F14394+5332), the wavelength of peak of the 
intensity (i.e. C80) is defined using the kinematic compo- 
nent at rest frame. In addition, in the case of the Syl IRAS 
F15130-1958, we used our modelling results to subtract the 
Fell emission prior to calculating the parameters. 

Figure [7| shows ^7 20 plotted against Kl and W2Q. 
INelson fc Whittle! (|l995l ) do not provide the center at 20% 
intensity level in their paper and therefore it is not possible 
to estimate the y4720 parameter for the objects in their sam- 
ple. At first sight, the figure shows the well-known tendency 
for the [OIII1 A5007 in Seyfer t gala x ies to have b l uewar d 
asymmetries (|Heckman et all Il98ll : IWhittle et ail . Il98&t ). 
In addition, the upper panel of the figure shows that, al- 
though the values of the different parameters found for the 
Sy-ULIRGs are consistent with the main body of points for 
the 4 comparison samples, a significant proportion (at least 
5 of 17, 30%) are outliers in these plots. Moreover, the lower 
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Figure 7. Upper panel: the asymmetry A/20 plotted against 
the kurtosis parameter Kl for the different sam ples of Sy galax- 
ies. Op en black square s: iHeckman et alj lll98ll). Open b lue dia - 
monds: I Whittle! (l985h . Open green triangles: IVeilleu 3 lll991al) . 
Filled red circles: our sample of Sy-ULIRGs. A/20 values higher 
(lower) than indicate blueward (redward) asymmetries, while 
Kl values lower (higher) than 1 indicate "peaky" ("stubby") pro- 
files, for a perfectly symmetric profile AI20 = 0, while Kl = 1 
for a Gaussian profile. Lower panel: the asymmetry A/20 plotted 
against W20. 



parison samples of Seyferts. Lower spectral resolution will 
generally re sult in larger lin e widths and lower asymmetry 
indices (see IVeilleuxl fl991bT ) . Since the A/20 values for the 
Sy-ULIRGs are higher on average than those for the com- 
parison samples, despite the possible resolution effects, the 
results in the table provide strong evidence that, in general, 
the emission lines of the Sy-ULIRGs are more asymmetric 
that those of non-ULIRGs Seyfert galaxies. Furthermore, 
the W20 values (and to a lesser extent the W50 values) for 
the Sy-ULIRGs are substantially higher than those of the 
objects in the comparison samples. Therefore, even allow- 
ing for resolution effectfl our results also suggest that the 
Sy-ULIRGs have broader emission lines compared to non- 
ULIRGs Seyfert galaxies. 

To further investigate the statistical significance of any 
differences between the kinematic properties of the Sy- 
ULIRGs and the non-ULIRG Seyferts we have performed a 
variety of stat istical te s ts. In the first instance, we used the 
2D KS test of H cacoc kl (|1983|) to investigate the significance 
of the differences between the positions of the Sy-ULIRGs 
and the comparison samples in Figure [7] The results are 
shown in Table [5] We find that, especially considering the 
comparison with all the Seyfert galaxy samples combined, 
the differences are significant at the > 3cr leveQ 

We have also compared the distributions of A/20, Kl, 
W50 and W20 separately using a ID two sample KS test. 
The results are consistent with those of the 2D KS test of 
iPeacockl <| 19831 ) and are shown in Table [5] We find that, al- 
though the emission line in Sy-ULIRGs are not particularly 
more "peaky" or "stubby" than those of the comparison 
samples, they tend to be more asymmetric, and they are 
certainly broader: the difference between the W20 and W50 
parameters for the Sy-ULIRGs in our sample and the ob- 
jects in the comparison samples is significant at the > 3<r 
level. 

At this stage is important to note that, while the com- 
parison samples include a significant fraction of broad line 
Seyfert 1 galaxies, our sample of Sy-ULIRGs includes only 
one such object. This could be a problem if the NLR kine- 
matics depend on Seyfert type (e.g. because of orientation 
effects). In this context, the results of lHeckman et all (|l98lf ) 
show that galaxies with a significant BLR component have, 
on average, narrower [OIII] lines that those that do not have 
such compone n ts. Ho weve r, the more compr e hensiv e stud- 
ies of I Whittle! (|l985l ) and iNelson fc Whittle! (|l995l ). found 
no evidence for differences between the line width distribu- 
tions of Syl and Sy2 galaxies of similar [OIII] luminosity 
and bulge properties. To explore any potential biases due 
to the different proportions of Syl objects in the Sy-ULIRG 
and comparison samples, we have repeated the analysis de- 
scribed above but this time including only in those objects 



panel in Figure [7] shows that Sy-ULIRGs have, in general, 
broader emission lines that the objects in the comparison 
samples. It is also apparent in the figure that all three of the 
objects with the most asymmetric profiles are Sy-ULIRGs. 

To investigate the possible differences between non- 
ULIRGs Seyfert galaxies and Sy-ULIRGs, Table[7]shows the 
mean values and estimated uncertainties for A/20, Kl, W20 
and VF50. We note that our observations of Sy-ULIRGs have 
a lower spectral resolution than those of most of the com- 



6 In reality, correcting for resolution has a negligible effect on 
the W20 measurements. For example, correcting a typical value 
W20 measured for the Sy-ULIRGs of W20 = 1000 km s" 1 for 
the 230 km s _1 instrumental width results in a decrease in W20 
of only 2.5%. 

7 Those objects included in more than one of the comparison 
samples are only considered once for this comparison. In partic- 
ular, we used the set of parameter values that are derived with 
the dataset for which the spectral resolution most resembles the 
spectral resolution of our spectroscopic data. 



22 J.Rodriguez Zaurin, C.N.Tadhunter 



Sample 


A/20 


Kl 


W20 

(km s _1 ) 


VK50 

(km s~ 1 ) 


Hcckman et al. (1981} 
N=36 


0.11±0.03 


0.88±0.02 


798±69 


455±40 


Whittle (1985) 
N=36 


0.10±0.02 


0.80±0.03 


632±40 


329±24 


Veilleux (1991a1 
N=16 


0.12±0.06 


0.78±0.04 


700±122 


374±73 


Nelson & Whittle (19951 
N=77 




0.87±0.02 


682±47 


388±28 


all samples 


O.lOiO.Ol 


0.85±0.02 


689±41 


384±16 


This paper 
N=17 


0.27±0.06 


0.74±0.06 


1423±131 


713±102 



Table 7. Mean and the standard error of the mean (<r/\/~N) for the kinematic parameters measured in the various samples. The entry 
in the table labelled as "all" correspond to the values considering the 4 comparison samples together. 



Sample 




A/20 vs Kl 


A/20 vs W20 


Sample 




A/20 


Kl 


W50 


W20 


(1) 


(2) 


(3) 


(4) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


Heckman et al. (19811 


D 


0.52 


0.56 


Heckman et al. (1981) 


D 


0.47 


0.43 


0.48 


0.63 




P(%) 


0.6 


0.2 




P(%) 


0.6 


2 


0.5 


<0.1 


Whittle (19851 


D 


0.46 


0.70 


Whittle (19851 


D 


0.53 


0.20 


0.62 


0.77 




P(%) 


2 


<0.1 


P(%) 


0.15 


65 


<0.1 


<0.1 


Veilleux (1991a1 


D 


0.36 


0.70 


Veilleux (1991a1 


D 


0.33 


0.18 


0.64 


0.75 




P(%) 


22 


<0.1 




P(%) 


24 


93 


0.1 


<0.1 


Nelson & Whittle (19951 


D 




0.62 


Nelson & Whittle (19951 


D 




0.33 


0.57 


0.66 




P(%) 




<0.1 




P(%) 




7.5 


<0.1 


<0.1 


All 


D 


0.43 


0.66 


all 


D 


0.47 


0.27 


0.56 


0.68 




P(%) 


0.6 


<0.1 




P(%) 


0.2 


17 


<0.1 


<0.1 



Table 8. The results of the 2D KS test of Peacock 1983, com- Table 9. The results of the ID two sample KS tests, comparing 

paring our sample of Sy-ULIRGs with the comparison samples. AI20, Kl, W50 and W20 for our sample of Sy-ULIRGs with 

Col 1: sample name. Col 2: 2D statistics. Col 3 and 4: the two tnosc Q f t h c four comparison samples, 
parameters considered for thc 2D KS test. 



classified as Sy2 for each individual sample. The results are 
consistent with the previous results described in this sec- 
tion, although with slightly less significance in the case of 
the A/20 parameter. 

To summarise, these results provide evidence that the 
emission line kinematics of the Sy-ULIRGs are more highly 
disturbed than those of the non-ULIRGs Seyferts in the 
comparison samples. Indeed, three of the objects in our sam- 
ple show more extreme emission line kinematics than any 
of Seyfert galaxies in the comparison samples. However, we 
emphasize that the size of our Sy-ULIRGs sample is rela- 
tively small. Therefore, the results from the K-S test must 
be interpreted carefully. Observations of larger samples of 
Sy-ULIRGs will be required to put this result on a firmer 
statistical footing. 



4-5.2 Comparison with the emission line kinematics of 
other types of AGN 

Whittle (1985) found evidence - albeit weak - for correla- 
tions between the emission line kinematics and the lumi- 
nosities of the associated of AGN. Given that the Seyfert 
galaxy comparison samples considered in the previous sec- 
tion contain a higher proportion of objects with low [OIII] 
luminosities than our Sy-ULIRG sample, it is therefore im- 
portant to make comparisons with samples than include 
higher luminosity objects. For this purpose we use the fol- 
lowing three AGN samples: the Palomar-Green (PG) sample 
of nearby quasars ( Boroson et al., 1982; Rose et al., 2012), 
the I Jin etT al. (2012) sample of unobscured Seyfert 1 nuclei 
and quasars, and the Rose et al. (2012) sample of nearby 
AGN selected from the 2MASS sample on the basis of their 
red near-IR colours. A major advantage of these samples 
compared to those used in the previous section is that the 
[OIII] emission line profiles have been fitted using the same 
multiple Gaussian fitting technique that we have used for 
the ULIRGs in this paper. The general properties of these 
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Sample Selection criteria Rcdshift range L[oiii] range Resolution References 

(W) (km s" 1 ) 

Sy-ULIRGs S 6 o >1 Jy, L IR > 10 12 L 0i z < 0.175 2 x 10 33 - 4 x 10 35 ~230 This paper 

emission line ratios 

PG Quasars UV excess, broad lines z < 0.5 5 X 10 33 - 5 X 10 36 ~350 Boroson et al. (1982') 

M B < -23 

Unobscured type I Type 1, X-ray bright z < 0.4 2x 10 33 - 3x 10 35 ~160 Jin et al. (20121 

Unobscured 

2MASS AGN (J-K) > 2 z < 0.29 3x 10 33 - 6x 10 35 ~550 Rose et al. 2012 

Table 10. Details of the properties of the higher luminosity AGN comparison samples. 



samples are compared with those of the Sy-ULIRGs in Table 

ESI 

We use the velocity shifts and FWHM of the broad and 
intermediate emission line components as our main indica- 
tors of the degree of kinematic disturbance. A complication 
with this approach is that individual objects — including 
many of the ULIRGs discussed in this paper — may have 
two or more broad and/or intermediate emission line compo- 
nents. In such cases we use the flux weighted mean shift and 
FWHM of the broad and intermediate emission line compo- 
nents (defined in Section 3.3); this yields single estimates of 
the shift and FWHM for each object and allows fair com- 
parison of the degree of kinematic disturbance between the 
various samples. Note that all the velocity shifts and FWHM 
are measured relative to the rest frames of the host galax- 
ies, and narrow emission line components (by our definition 
FWHM < 500 km s _1 ) are not included in the analysis. 

Figure [8] shows the velocity shift plotted against the 
FWHM for the Sy-ULIRGs and the three comparison sam- 
ples. Again, the figure shows the well known tendency for the 
kinematically disturbed components to be blueshifted rela- 
tive to the host galaxy rest-frames. Also, there is a tendency 
for the objects wit h the largest blueshif ts to have the largest 
FWHM (see also ISpoon fc HoltJ . [20091 ') , although the con- 
verse in not true, indeed the object with the largest FWHM 
- the Sy-ULIRG F23389+0300 — shows a slight redshift. 
Consistent with the results of the previous section, we find 
that, although a significant proportion of the Sy-ULIRGs 
fall within the main envelope of points defined by the three 
comparison samples (AV > -500 km s"\ FWHM < 1500 
km s _1 ), 50% are outliers on these plots with AV <-500 
km s _1 and/or FWHM > 1500 km s" 1 ; this proportion is 
higher that found in the PG quasar (15%), unobscured type 
I (12%), and 2MASS (34%) samples. We further note that 
detailed examinations of the spectra of the objects in the 
high luminosity comparison samples reveals no objects with 
the type of highly structured, extremely blueshifted, multi- 
ple component [OIII] emission profiles that are observed in 
the Sy-ULIRGs F134 51+1232 and F 14394+4332. 

We again use the lPeacockl (|l983h 2D KS test to investi- 
gate the significance of the differences between the positions 
of the Sy-ULIRGs and the comparison samples in Figure 
[8] The results are shown in Table 1111 From this it is clear 
that we can reject the null hypothesis that the Sy-ULIRGs, 
PG quasars and unobscured type I AGN objects are drawn 
from the same parent population at the ~2% (~ 2cr) level 
of significance, but do not find a statistically significant dif- 
ference between the Sy-ULIRGs and the local 2MASS AGN 
population. 



Comparison sample N 2D KS statistic (D) P(%) 

PG Quasars 85 0.45 1.9 

Unobscured type I 30 0.47 2.4 

2MASS AGN 30 0.32 38 

Table 11. The results of the 2D KS test of Peacock 1983, compar- 
ing our sample of Sy-ULIRGs with the three comparison samples. 



FWHM Shift 



Comparison sample 


D 


P(%) 


D 


P(%) 


PG Quasars 


0.48 


0.3 


0.34 


7.6 


Unobscured type I 


0.43 


1.9 


0.40 


4 


2MASS AGN 


0.28 


38 


0.24 


(50 



Table 12. The results of the ID two sample KS tests, comparing 
the FWHM and velocity shift distributions for our sample of Sy- 
ULIRGs with those of the three comparison samples. 



We have also compared the distributions of FWHM and 
velocity shift separately using a ID two sample KS test. The 
results shown in Table E2 reveal that the most significant 
differences are between the distributions of FWHM for the 
Sy-ULIRG and the PG Quasars, with the difference signif- 
icant at the 0.3% (~ 3a) level. This difference is consistent 
with that found between the line widths of the Sy-ULIRG 
and non-ULIRG Seyfert comparison samples in the previ- 
ous sectior[f|. The differences between the ID distributions 
of velocity shift for the Sy-ULIRGs, PG quasars and Un- 
obscured Type I samples are less significant and, consistent 
with the results of the 2D KS test, the ID tests provide 
no evidence for significant differences between the FWHM 
and velocity shifts when comparing the Sy-ULIRGs with the 
local 2MASS sample. 

Together, the results of this section show a tendency for 
the Sy-ULIRGs to display a higher degree of kinematic dis- 
turbance than the PG quasars and unobscured type I AGN, 
but again observations of larger samples of Sy-ULIRGs will 
be required to strengthen the statistics. 



8 Note that, in this case, there is no issue with differences between 
the spectral resolutions of the samples, because the FWHM values 
have been corrected for instrumental resolution. 
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Figure 8. Flux weighted mean shifts plotted against the flux 
weighted mean FWHM for the ULIR Gs in our sample ( black-filled 
circles) and the PG QSOs (top), the |jin et al,l d2012h sample of 
unobscured Seyfert 1 nuclei and quasars (middle) and the Rose 
et al. (2012) sample of nearby AGN selected from the 2MASS 
sample on the basis of their red near-IR colours (bottom). The 
figure show that the proportion of objects that show extreme 
emission line kinematics in the sense that they are outliers in the 
plot is higher for the Sy-ULIRGs than it is from the comparison 



5 CONCLUSIONS AND FUTURE WORK 

Most of the recent investigations of outflows in ULIRGs 
have focussed on observations of the neutral and molec- 

n 



Fischer et al.l . |2010| ; 
20 111 ). Thesehave re- 



ular gas fe.g | Rupke et al.l. l2005ct 



1 r 



IRupke fc VeilleuxtboTll ; ISturm et al 
vealed massive, highly energetic outflows in several objects, 
supporting the view that the feedback effect associated with 
the outflows can have a significant impact on the evolution 
of the host galaxies. However, there remain questions about 
whether these molecular and neutral outflows are driven by 
the circum-nuclear starbursts or by the AGN that are found 
in the nuclei of some of the sources. Indeed many such stud- 
ies have failed to find significant differences between the 
outflows in ULIRG s with and without powerful AGN (e.g. 
IRupke et ail . l2005d ) . 

Given that (a) we find that a significantly higher frac- 
tion of the ULIRGs with optical Seyfert nuclei show evi- 
dence for nuclear outflows than those lacking such nuclei, 
and (b) the outflow components have line ratios consis- 
tent with AGN- rather than starburst-photoionization, it 
is highly likely that the nuclear outflows in the Sy-ULIRGs 
are indeed driven by the AGN. We also find that the Sy- 
ULIRGs present outflows that are are more extreme than 
those in samples of local non-ULIRG Seyfert galaxies. 

Clearly there is now evidence that the AGN associated 
with ULIRGs drive warm outflows. What remains uncertain 
is whether these outflows are sufficiently energetic and mas- 
sive to affect the star formation histories and evolution of 
the host merger remnants. In section 3.8 we estimated the 
ranges of mass outflow rates and kinetic powers that are con- 
sistent with our observations. At the upper end of the range, 
the mass outflow rates and kinetic powers of the warm out- 
flows in Sy-ULIRGs are comparable with those typically es- 
timated for the neutral outflows detected using the optical 
NalD li nes (6 < M < 400 M pyr" 1 ; 3 x 10 40 < E < 10 44 
erg s-^ lRupke et al.l . l2005bl lc]; lRupke fc Veilleuxl, 120111 ) and 
the molecula r outfl ows using the hydroxyl molecule OH 
l|Sturm et all 1201 ll . up to 1.2xl0 3 M yr" 1 ), and would 
be likely to have a major impact on the host galaxies. On 
the other hand, if the lower limiting estimates were to prove 
more typical of the true outflow properties, the warm out- 
flows would be insignificant compared with the neutral and 
molecular outflows, and therefore unlikely to significantly 
affect the evolution of the host galaxies. 

Despite the high degree of uncertainty, there is cause 
for optimism that the mass outflow rates and kinetic power 
estimates will be substantially improved in the near future. 
In particular, using the transauroral [Oil] and [SII] line ra- 
tios it will be possible to make density and re ddening es- 
timat es that are considerably more precise (see lHolt et al.l . 
2011); the outflow radius estimates will benefit from HST 
narrow-band (e.g. iBatcheldor et all . 120071 ) and AO-assisted 
near-IR imaging; and the reddening (and hence H/3 lumi- 
nosities) will be further refined by combining optical and 
near-IR measurements of the hydrogen recombination line 
fluxes. 
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APPENDIX: Description of individual objects 

F00188-0856. Veilleux et aD (|l999h classified this compact 
object as a LINER a t optic al wavelengths. However, based 
on the iKewlev et al.l (|2006l ) (hereafter Kwl06) classification 
scheme, Yuan et al. (2010) classified it as "Sy2:", meaning 
that F00188-0856 has a Sy2 spectral type in 2 of the 3 BPT 
diagrams use for its optical classification (see Yuan et al., 
2010 for details). Therefore, we included this object in the 
sample discussed in this paper. 

F00188-0856 is one of the two cases (along with 
F23327+2913), for which the [OIII] AA4959, 5007 and H/3 
emission lines have a low equivalent width, and are poten- 
tially affected by noise and residual structure in the under- 
lying continuum. Therefore these lines are not, in this case, 
suitable for deriving a model to fit to the other lines in the 
spectrum. For this object, we used the approach of first fit- 
ting the Ha+[NII] complex and then modelling the other 
emission lines in the spectrum with the same components 
and the same velocity widths and shifts as Ha. In addition, 
due to the presence of a telluric absorption feature that co- 
incides with the location of the [SII]AA6716,6731, it was not 
possible to model this emission blenrfl The results of the 
Ha+[NII] modelling are shown in Figure [T] 

The best fitting model comprises two components: a 
narrow component with FWHM = 269 ± 5 km s _1 plus 
an intermediate component with FWHM = 904 ± 75 km 
s _1 blueshifted by -253 ± 36 s" 1 with respect to the nar- 
row component. Using these results is possible to measure 
the line ratios associated with each individual component 
and then study the corresponding ionization mechanisms us- 
ing the already mentioned Kwl06 classification scheme. We 
find that the line ratios associated with the narrow compo- 
nent are consistent with a composite spectral type (CP) or 
LINER, while those corresponding to the intermediate com- 
ponent are consistent with Sy2 ionization, but lie close to 
the Sy2/LINER limiting region. 

F01004-2237. The optical spectrum of this single nu- 
cleus system shows a mixture of HII and Sy2 features 
I Veilleux et all 1 19991 ). However, it was classified as a "Sy2:" 
in the work of Yuan et al. (1010) and therefore, included it 
in our sample. 

Figure [1] shows the [OIII] model for this galaxy, which 
includes three components: an unresolved, narrow compo- 
nent, an intermediate component (FWHM = 850 ± 71 km 
s _1 ) blueshifted by -230 ± 32 km s _1 with respect to the 
narrow component and a broad (FWHM = 1590 ± 92 km 
s _1 ) component blueshifted by -999 ± 80 km s _1 with re- 
spect to the narrow component. Note that the broad, shifted 
components dominate the [OIII] emission line flux in this 
case. 

Due to the number of components, the degeneracy be- 
tween the different models that adequately fit the emission 
line profiles is a relatively important issue for this object. 



9 There are cases for which a telluric absorption feature af- 
fects the profiles of some lines (mainly [SII] AA6716,6731 and the 
Ha+[NII] AA6548,6583 complex) but it is still possible to model 
the emission lines. However, the flux values derived from the mod- 
elling will be relatively unconstrained. An example of the typ- 
ical effect of a telluric absorption feature on the profile of the 
Ha+[NII]AA6548,6583 complex in shown in Figure ITU 



For example, when the intensities of the different compo- 
nents were left as a free parameter during the modelling of 
the Ha+[NII]complex, the [Nil] A6583/Ha line ratio corre- 
sponding to the broad component was ~4. As shown in Yuan 
et al. (2010), their Figure 2, this line ratio is always smaller 
than ~3. Therefore, during the modelling, we set the value 
of the [NII]A6583/Ha line ratio to 3, obtaining a good fit to 
the [NIi]+Ha blend. 

In addition, the [SII]6716/6731 line ratio was higher 
than 1.42 (the low density limit) for the narrow component, 
but lower than 0.44 (the high density limit) for the two broad 
components (i.e., the results were unphysical) when making 
free fits to the [SII] blend. Therefore, during the modelling, 
we also forced the values of these line ratios to be 1.42 and 
0.44 for the narrow and the two broad components respec- 
tively. Then, the value of one of these line ratios was changed 
in steps on 0.02 while the values of the other two remained 
unchanged until the fit was deemed not adequate. We es- 
timate a lower limit for the intensity ratio [SII]6716/6731 
of 1.2 for the narrow component while acceptable fits are 
obtained for any line ratio value within the allowed range 
(0.44 < [SIIJ6716/6731 < 1.42) in the cases of intermediate 
and the broad components. 

Regarding to the ionization mechanisms associated with 
each of the individual components, we find that the line ra- 
tios corresponding to the narrow and intermediate compo- 
nents are consistent with HII ionization or fall close to the 
HII/Sy2 limiting region in the diagrams. On the other hand, 
the line ratios associated with the broad component are con- 
sistent with NLAGN and Sy2 ionization in two of the BPT 
diagrams, while they fall close to the HII/Sy2 limiting re- 
gion in the [SII] diagram. Given these results, it is clear why 
the spectral classification of this object ba sed on line ra- 
tios d erived from integrated line fluxes (e.g. IVeilleux et al.l . 
1999), was ambiguous in previous studies; our analysis shows 
the clear presence of an AGN associated with the broad, 
blueshifted emission line components. Note that the ap- 
proach used to model the [SII]6716/6731 doublet does not 
significantly change these results. 

The nuclear emission line kine matics of this objec t 
have been d i scusse d in t he past by | Armus et all |l988), 
iLfpari et all (l2003h. and IFarrah et all d2005l) . Whi le the 
spectra of Armus et alj (|l988l ) and lLfpari et alj ((2003) show 
the same blue wings to the [OIII] emission line s that we de- 
tect in the spectrum presented in this study, IFarrah et al.l 
(2005) report that they fail detect the blue wings in their 
HST/STIS spectrum (0.2" slit width) that isolates the com- 
pact knot in the nucleus of the galaxy. This is interesting be- 
cause it suggests that the AGN outflow may be significantly 
extended in this object. However, the substa ntially worse 
spectr al resolution of the HST/STIS spectra of IFarrah et alj 
(2005) compared with the spectra presented in this paper 
may have prevented the detection of the shifted features. 

F12072-0444. As in the case for F00188-0856, the pres- 
ence of a telluric absorption feature coinciding with the lo- 
cation of the [SII]AA6716,6731 blend prevents any attempt 
to model this feature for this Sy2 galaxy. The best fitting 
model to all the other observed emission lines comprises 3 
components: a narrow component (FWHM = 275 ± 8 km 
s _1 ) and intermediate component with FWHM of 525 ± 28 
and blueshifted -276 ± 30 km s _1 with respect to the narrow 
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component and a broad component with FWHM of 1343 ± 
81 km s" 1 and blueshifted -446 ± 35 km s" 1 . The line ra- 
tios measured for all of the three components are consistent 
with NLAGN or Sy2 ionization. 

F13305-1739 F13305-1739 has by far the highest emis- 
sion line luminosity of all the objects in our sample (see 
Table 1); its [OIII] emission line luminosity would lead to 
it being classified as a type II q uasar according the the cri- 
teria of IZakamska et all (|2003l ) . and it shows a rich emis- 
sion line spectrum from the UV ([MgII]A2798) to the red 
end ([SII]AA6716,6731) of its spectrum. Figure [T] shows the 
[OIII] model that successfully fits all the other emission lines 
detected in the optical spectrum of this source. Clearly, the 
broad/shifted components dominate the [OIII] emission line 
flux in this case. 

We use the case of F13305-1739 as an example of our 
modelling technique, and Figure [5] shows the modelling re- 
sults obtained for all the emission lines using the so called 
[OIII] model. This model comprises three comp: a narrow 
component (FWHM=435±25 km s _1 ), a first broad com- 
ponent (Bl) with FWHM of 1276±29 km s _1 blueshifted 
by only -36 ± 14 km s _1 with respect to the narrow compo- 
nent and a second, broad component (B2) with FWHM of 
1685±286 km s" 1 and blueshifted by -281±29 km s" 1 with 
respect to the narrow component. Figure [U] shows that, in 
general, the Bl component dominates the flux contribution 
for the majority of the emission lines detected. 

As in the case of F01004-2237, degeneracy is an impor- 
tant issue in the case of F13305-1739. When the intensities 
of the different components were left as a free parameter dur- 
ing the modelling the [SIIJ6716/6731 line ratio was higher 
than 1.42 (the low density limit) for the narrow component, 
but but this ratio was lower than 0.44 (the high density 
limit) for the two broad components (i.e., the results were 
unphysical). Using the same approach described before, we 
estimate a lower limit for the intensity ratio [SII]6716/6731 
of 1.2 for the narrow, component and an upper limit of 0.6 
for the broad components. 

In terms of the ionization mechanisms, the line ratios 
measured for all three kinematic components are consistent 
with Sy2 ionization. Note that the approach used when mod- 
elling the [SII]AA6716,6731 doublet does not significantly 
change these result. 

F13428+5608 (Mrk273). This well-known ULIRG is the 
closest in our sample of Sy-ULIRGs (z — 0.037), and in 
optical images shows an impressive tidal tail that extends 
over 30 kpc towards the south of the g alaxy. A dou b le nu- 
cleus structure is vis i ble at near-IR jArmus et al.l. I199C ; 



Maiewski et all . Il993l; iKnapen et al.l , ll997 



Scoville et al 
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Cole et al.l , ll999T l wavelengths, with a nuclear separation of 
~ 700 pc, although the two nuclei are unresolved in our 
spectroscopic data. 

As shown in Figure [T] the best fitting model for this 
galaxy includes two components: a narrow component with 
a FWHM of 449±7 km s" 1 plus a broad component with a 
FWHM of 1369±62 km s _1 that is redshifted by 14±12 km 
s _1 with respect to the narrow component. 

We find that, when the intensities of the two compo- 
nents are left as free parameters during the modelling of the 
Ha+[NII] complex, the [NIL] A6583/Ha line ratio measured 



for the broad component (FWHM = 1369 km s _1 ) is signif- 
icantly higher than 3, which, based on the results of Yuan 
et al. (2010), can be considered as an upper limit for this 
line ratio. In addition, the [SII]6716/6731 line ratio obtained 
for the broad component is outside the allowed range (i.e. 
[SII]6716/6731 < 0.44, the high density limit). Therefore, as 
discussed before, we set the value of the [Nil] A6583/Ha line 
ratio to 3 and the [SII]6716/6731 line ratio at the high den- 
sity limit. Then we gradually increased the value of the 
[SII]6716/6731 line ratio until the overall fit was poor, which 
occurred at the value of 0.8. 

In terms of the ionization of the gas, the line ratios 
measured for the narrow component lie close to the limit- 
ing regions between CP and narrow line AGN (NLAGN) 
ionization, and LINER and Seyfert2 galaxies in the Kwl06 
classification scheme. On the other hand, the line ratios mea- 
sured for the broad component are consistent with Sy2 ion- 
ization. Note that the approach used when modelling the 
Hq+[NII] complex and the [SII]AA6716,6731 doublet does 
not change this result. 

F13451+1232. A detailed discussion of the spectacu- 
lar emission line kin ematics thi s obje ct can be found in 
iHolt et al l (|2003l ) and lHolt et ail |201ll ). 

F14394+5332. F14394+5332 shows the most spectacular 
[OIII] emission line kinematics of all the objects considered 
in this paper (see Figure 1), with the clear detection of three 
distinct narrow components, as well as an underlying broad 
component. Based on our analysis of the spatially extended 
[OIII] emission line kinematics in this object, the redmost 
narrow component in the integrated spectrum represents the 
galaxy rest frame. In this case, the two other narrow com- 
ponents are blueshifted by -700 and -1360 km s _1 relative 
to the rest frame, and the underlying broad components are 
also blueshifted; clearly this is a case in which the compe- 
nents that are broad and/or blueshifted dominate the [OIII] 
emission line flux. The multi-componen t nature of the [QU I] 
lines in this object was also detected bv lLfpari etakl l|2005t) . 
albeit at lower S/N. 

The profiles of the [OIII] emission lines in this object 
are strikingly different from those of most of the other emis- 
sion lines. This is illustrated by the comparison between the 
H/3, [OIII] and [NI] emission lines shown in Figure 1101 In 
particular, we do not clearly detect the blueshifted narrow 
lines in any of the lower ionization lines such as [OII]A3727, 
H/3, [NLU5200, [OI]AA6300,636£3, Ha, [NII]AA6548,6584 
and [SII] AA6717,6731, apart from possible hints of the first 
shifted narrow component in the form of a "shelf" in the 
blue wings of [NI]A5200 and [OI]A6300. On the other hand, 
the [NeIII]A3869 line has a profile that is consistent with 
that of [OIII]A5007, although the lower S/N, lower effective 
velocity resolution, and possible residual features from the 
subtraction of the (complex) underlying stellar continuum, 
preclude an accurate comparison between the line profiles. 
We further note that the [NeV]A3425 and HeII(4686) fea- 
tures show evidence for broad, blueshifted features, but the 
S/N is again too low to allow a detailed comparison with 
the [OIII] line profile. 



10 N.B. The blue wing of the [OI]A6300 line may be significantly 
affected by telluric atmospheric B-band absorption. 
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Figure 9. The results obtained by using the [OIII] model to fit all the optical emission lines detected in the spectrum of Ff3305-f 739. 
The resulting model is overplottcd (bold line) on the extracted spectrum (faint line) . The different kinematic components are also plot in 
the figure. The solid, dashed and dotted line correspond to the narrow component and the two broad components respectively (referred 
as N, Bf and B2 in the text). The figure shows that Bl (dashed line) dominated the flux contribution for the majority of the emission 
lines detected in the optical spectrum of the source. 



The fact that the discrete blueshifted [OIII] emission 
lines components are not clearly detected in the other lines, 
many of which have similar total emission line fluxes, im- 
mediately implies that the gas emitting these blueshifted 
narrow components has a high ionization state. Nonetheless, 
the low ionization lines do show some evidence for kinematic 
disturbance in the form of blue wings. Overall, we find that 
we can obtain an acceptable fit to all the low ionization and 
Balmer lines in the nuclear spectrum using a model that in- 
cludes a narrow component with the same kinematic prop- 
erties as those of the redmost (rest frame) narrow compo- 



nent detected in the [OIII] line, plus two broad components 
blueshifted shifted by -156 and -1574 km s" 1 with FWHM 
of 1272 and 1401 km s _1 respectively. We also find that we 
can fit the [OIII] profiles with the same two broad compo- 
nents plus the three narrow components described above; 
this is fit is shown in Figure 1, and all the corresponding fit 
components are listed in Table 2. 

To get some idea of the ionization of the kinematically 
disturbed emission line components, we have compared the 
sum of the fluxes of all the blueshifted components of [OIII] 
(2 narrow + 2 broad) with the sum of the fluxes of the 
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Figure 10. The modelling results obtained for [OIII], H/3 and [NI]A5200. The resulting model is overplotted (bold line) on the extracted 
spectrum (faint line). The different components are also plot in the figure. The solid, dashed and dotted line correspond to the narrow 
components (3 in the case of the [OIII] emission line referred as Nl, N2 and N3 in Table 2) and the two broad components respectively 
(referred as Bl and B2 in Table 2). The figure shows that the blucshiftcd narrow lines (N2 and N3) are not clearly detected in the lower 
ionization lines such as H/3 or [NIJA5200. 



blueshifted components fitted to the lower ionization lines 
(2 broad). The resulting line ratios are labelled "B" in Table 
3. Interestingly, despite the fact that the discrete blueshifted 
narrow components in [OIII] must be emitted by gas with a 
very high ionization state, the overall B line ratios are low 
ionization in character, with [Oil]/ [OIII] ~ 1 and [OIII] /Ha 
~ 1, although more consistent with AGN than starburst 
heatin j"! On the other hand, the line ratios for the red 
narrow component fall squarely in the HII/AGN region of 
all three diagnostic diagrams. 

Clearly, the extraordinary kinematical and ionization 
properties of this source warrant further investigation with 
higher resolution and S/N optical spectroscopy. 

F15130— 1958. Three components were required to ade- 
quately model the profile of the emission lines observed in 
the optical spectrum of this single nucleus ULIRG. A first, 
intermediate component (II, FWHM = 545±69 km s -1 ) 
component, a second, intermediate component (12, FWHM 
= 700±188 km s" 1 ) blueshifted by -351±280 km s _1 with 
respect to II, and a broad (FWHM = 1630±41 km s^ 1 ) 
component blueshifted by -725±131 km s _1 with respect to 
the narrow component. 

As in previous cases, when only the shifts between the 
different components and their widths of are fixed during 
the modelling of the Ha+[NII] complex, the results obtained 
lead to unphysical Ha/H/3, [NII]A6583/Hq line ratio values 
for 12, and [SII]6716/6731 line ratios for 12 and the broad 
component. Therefore, during the modelling, the H/3 flux of 
the 12 component was set such the the Ha/H/3 line ratio as- 
sociated with this component is equal or higher than that of 
II (Ha/H/3 = 3.74 and 4.01 for II and 12 for the best fitting 
model). In addition, the [NII]A6583/Hq was set to 3 and the 
[SII]6716/6731 line ratio was set to 1.42, 0.44 and 0.44 for 
II, 12 and the broad component respectively. The values of 
each of the [SIIJ6716/6731 line ratios was then changed in 
steps on 0.02, while the values of the other two remained un- 
changed, until the fit was deemed poor. Hence, we estimate 
a lower limit for the intensity ratio [SIIJ6716/6731 of 1.2 for 



Note that integrated H/3 flux for the kinematically disturbed 
components is likely to be strongly affected by the underlying 
H/3 absorption. Therefore B ratios measured relative to H/3 are 
unlikely to be reliable 



II, while acceptable fits are obtained for any line ratio value 
within the allowed range (0.44 < [SII]6716/6731 < 1.42) in 
the cases of 12 and the broad component. 

Finally, we note that the line ratios of II are consis- 
tent with NLAGN ionization or fall close to the LINER/Sy2 
limiting region in the optical diagrams, whereas the values 
measured for the 12 and broad components are consistent 
with Sy2 ionization. 

F15462-0450. This is the only Seyfert 1 galaxy included 
in the sample discussed in this paper. Figure 1 shows the 
best fitting model obtain from fitting simultaneously the H/3, 
[OIII]AA4959,5007 and the broad FeIIAA4924,50l£3 emis- 
sion lines, typical of Syl galaxies. Due to the large number 
of emission lines close together in wavelength and bearing in 
mind that multiple (and broad) kinematic components con- 
tribute to each individual emission line, clearly overcoming 
the degeneracy issue is not easy in the case of F15462-0450. 

Overall, we find that the Balmer lines (H<5, H7,H/3 and 
Ha) are adequately fitted using three kinematic components: 
a narrow component of width 147 ± 47 km s -1 at rest 
frame, plus a first broad component of width 1509 ± 47 km 
s _1 redshifted 169 ± 8 km s _1 with respect to the narrow 
component and a second broader component of width 3673 
± 137 km s" 1 redshifted 292 ± 21 with respect to the 
narrow component. 

On the other hand, the forbidden lines 
([NeV]AA3346,3425, [OII]AA3726,3728, [NeIII]AA3868,3968, 
[OIII]AA4959,5007, [NII]AA6548,6583 and [SII]AA6716,6731) 
are adequately fitted using two kinematic components that 
correspond to those presented in Table 2, and labelled as 
the [OIII] model for this galaxy: a narrow component at 
rest frame, consistent with that found for the Balmer lines 
(i.e. FWHM = 147 ± 47 km s -1 and referred as 'N' in 
Table 2) plus a broad component (FWHM = 1426 ± 47 
km s _1 ) blueshifte by -822 ± 25 km s _1 with respect to 
the narrow component (referred to as 'B' in Table 2). Note 
that the broad component dominates the [OIII] emission 
line flux in this object. 

Note that the broad component associated to the for- 
bidden lines (B) it is still required to adequately model the 

12 Note that, in order to model the broad Fell lines, we assumed 
that they have the same profile as the broad H/3 line. 
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Figure 11. In the case of F16156+1046 it was not possible to 
model all the emission lines using a single model. However, with 
the exception of the [OIII] AA4959,5007, it was possible to model 
the rest of the observed emission lines with the model shown in 
the figure. This model comprises an unresolved component at rest 
frame (solid line), a narrow component (FWHM = 446 ± 86 km 
s _1 ) rcdshifted 33 ± 29 km s" 1 with respect to the unresolved 
component (dashed line) and a third, broad component (FWHM 
= 1225 ± 108 km s" 1 ) that is blueshifted by -150 ± 48 km 
s — 1 relative to the narrow, unresolved component (dotted line). 
The resulting model (bold line) is ovcrplotted on the extracted 
spectrum. 

profile of the H/3 emission line. Such component is repre- 
sented by the most redshifted, dashed- line Gaussian in the 
H/3 fit in Figure 1. Most likely, this component is also present 
in the other Balmer lines. However, the dominant flux con- 
tribution associated to the emission from the Broad Line 
Region (BLR) as well as the relatively low S/N in the cases 
of the H<5, H7 emission lines prevents a clear detection of 
this kinematic component for the other Balmer lines. In ad- 
dition, due to the dominant BLR emission, it is not possible 
to constrain the flux contribution of the broad component 
to the [NII]AA6548,6583 emission lines with any accuracy. 

Given the fact that the broad component present in 
the forbidden lines is not clearly detected in the Ha emis- 
sion line, no attempt to study the ionization mechanisms 
was made for that component. Interestingly, the line ratios 
measured for the narrow component detected in all lines are 
consistent with HII ionization. 

F16156+0146. For the study presented here we 
used the 5kpc-I aper t ure d escribed in the work of 
iRodrfguez Zaurm et al.l (|2009h . This aperture samples 
the north-western nucleus of this double nucleus system, 
which is associated with the AGN activity. This is one of 
the four objects in our sample (along with F13451+1232, 
F14394+5332 and F23389+0300) for which is not possible 
to find one model that fits all the emission lines observed 
in the optical spectrum. However, it is possible to model 
all the emission lines using two different approaches. A 
first approach in which a different model is used for each 
of the observed optical spectrum, and a second one in 
which we use two models: one that adequately fits the 
[OIII]AA4959,5007 emission lines, and a second model that 
fits all the other optical emission lines well. 

The best fitting model for [OIII] AA4959,5007, which is 
the same for the two approaches, is shown in Figure [1] and 
comprises three components: an unresolved narrow compo- 
nent, an intermediate component (FWHM = 804 ± 20 km 
s _1 ) blueshifted -186 ± 10 km s _1 with respect to the nar- 
row component and a broad (FWHM = 1535 ± 59 km s _1 



) blueshifted by -374 ± 26 km s _1 with respect to the nar- 
row component. The second model, shown in Figure [TT1 for 
the H/3 line comprises another 3 components: an unresolved 
component, a narrow component (FWHM = 446 ± 86 km 
s _1 ) redshifted by 33 ± 29 km s _1 with respect to the unre- 
solved component, and a third, broad component (FWHM 
= 1225 ± 108 km s _1 ) that is blueshifted by -150 ± 48 km 
s _1 relative to the narrow, unresolved component. 

Overall, the fact that we do not find one adequate 
model for all emission lines suggest that, instead of a discrete 
number of kinematic components, there may be a continu- 
ous range of density, ionization and kinematics through the 
Narrow Line Region (NLR), leading to each line having a 
different profile. We did not attempt to estimate line ra- 
tios associated with the different components in the case of 
F16156+0146. 

F17044+6720 . This ob j ect w as originally classified as a 
LINER bv lVeilleux et alj i| 19991 ). However, optical spectro- 
scopic classification of this source was recently reviewed 
by Yuan et al. (2010), who found that the line ratios of 
F17044+6720 are typical of a Sy2 galaxy. 

Figure [1] shows the [OIII]AA4959,5007 model for this 
galaxy comprising two components: a narrow component of 
width FWHM = 290 ± 5 km s _1 plus a broad component 
of width FWHM = 1765 ± 100 km s" 1 that is blueshifted 
by -553 ± 65 km s _1 with respect to the narrow component. 
Due to the presence of a telluric absorption feature it was 
not possible to model the [SII] AA6716,6731 emission lines for 
this galaxy. The line ratio values measured for the narrow 
component are consistent with CP ionization, or lie close to 
the limiting region between LINER and Sy2, while the line 
ratios associated with the broad component are consistent 
with NLAGN, or fall close to the LINER/Sy2 ionization. 

F171 79+5444. Three components are required to model 
the optical spectrum of this single nucleus Sy2 galaxy. Figure 
[T]shows an [OIII] model including a narrow component with 
a FWHM of 356 ± 22 km s _1 , an intermediate component 
with a FWHM of 515 ±51 km s _1 that is redshifted by 
246 ± 67 km s _1 with respect to the narrow component, 
and a third, broad component with a FWHM of 1562 ± 
41 km s _1 that is blueshifted with respect to the narrow 
component by -119 ± 61 km s _1 . The line ratios measured 
for the narrow component are consistent with CP or LINER 
ionization, while the line ratios of the shifted components 
are, in general, consistent with Sy2 ionization. 

However, for this galaxy, it is important to add a caveat 
about the uniqueness of the best fitting model. As described 
above, there are two close narrow components that have rel- 
atively similar widths (356 ± 22 km s _1 and 515 ± 51 km 
s" 1 ). Therefore, it is possible that this narrow line splitting 
represents unresolved rotation, or traces a large-sca l e bipo - 
lar outflow in the gas (see, for example iHolt et all [2008). 
To test this scenario we set the widths of the two narrow 
components to have the same value during the modelling of 
the [OIII] AA4959, 5007 emission lines; the result is presented 
in Figure 1121 The model, that also successfully fits all the 
observed emission lines, comprises three components: a nar- 
row component (Nl, FWHM = 386 ± 14 km s" 1 ), a second 
narrow component (N2) with the same width redshifted 301 
± 21 km s" 1 with respect to Nl and a broad component 
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Figure 13. Fits to the optical emission lines in F23060+0505. Top panel: H/3and Ha+[NII] AA6549,6583 fit including the two narrow 
components (Nl and N2) plus the broad components corresponding to the [OIII] model for this galaxy. Lower panel: H/3 fit including 
the already mentioned three components plus an underlying broad component (representing the BLR). When only 3 components are 
included in the modelling, the blue, and paricularly the red wings of the line profiles are not adequately fitted. The problem disappears 
when including the underlying broad component. 



(FWHM = 1543 ± 34 km s _1 ) blueshifted by -157 ± 33 km 
s _1 with respect to Nl. 

Note that when only the widths and shifts of the dif- 
ferent components were constrained for the two best fitting 
models described before, the Ha/H/3 line ratios associated 
with the intermediate and the N2 component were lower 
than 2.88, the theoretical lower limit assuming case B re- 
combination theory in absence of reddening. A similar is- 
sue is found for the broad component in the case of the 
[SII]AA6716,6731 doublet, where the [SII]6716/6731 line ra- 
tio is found to be higher than 1.42, the low density limit. 
This is due to the fact that the flux contribution in the case 
of the [SII]AA6716,6731 emission lines is dominated by the 
narrow components at rest frame (for the two best fitting 
models considered). Therefore, the intensities of the broad 
components are relatively unconstrained. 

To overcome these problems the H/3 flux contribution 
of the intermediate and the N2 was forced to have an in- 
tensity such as the Ha/H/3 line ratio associated with this 
component is higher or equal to that of the narrow compo- 
nent at rest frame. In addition, following the approach used 
for other objects, the [SII]6716/6731 line ratio of the shifted 
components were set to 0.44 during the modelling and then 
slightly increased in steps of 0.02 until the fit became poor. 
As we mentioned before, the flux contribution of these com- 



ponents is relatively small, and adequate fits are obtained 
for any value of this ratio within the allowed range (0.44 
>[SII]6716/6731>1.42). 

In terms of line ratios we find similar results for the two 
models described above. The values measured for the nar- 
row component is consistent consistent with CP and LINER 
ionization. On the other hand, the line ratios of the I and 
B components are consistent with Sy2 ionization. We note 
that telluric atmospheric absorption features affect, to some 
extent, the profile of the [OIJA6300 emission line and the 
blue wing of the [SII]AA6716,6731 doublet. Therefore, the 
line ratios for the broad components of the latter lines are 
relatively unconstrained for either of the models considered 
here. 

F23060+0505. Strong emission lines are observed from the 
UV end (MgIIA2798) to the red end ([SII]AA6716,6731) of 
the spectrum of this compact galaxy. The best fitting model 
presented in Figure[T]comprises 4 components: a narrow (Nl, 
FWHM = 376 ± 6 km s~ ), a second, unresolved, narrow 
component (N2) redshifted by 306 ± 4 km s" 1 with respect 
to the narrow component at rest frame, a broad component 
(FWHM = 1001 ± 23 km s" 1 ) blueshifted by -310 ± 14 
km s _1 with respect to the narrow component at rest frame 
and a very broad component (FWHM = 2148 ± 125 km s" 1 
) that is blueshifted -1073 ± 122 km s _1 with respect to 
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Figure 12. A second, possible [OIII] model that adequately fits 
all the observed emission lines in the case of F17179+5444. The 
different components shown in figure are: a narrow component 
(Nl, FWHM = 386 ± 14 km s -1 ), a second narrow component 
(N2) with the same width redshifted 301 ± 21 km s — 1 with re- 
spect to Nl and a broad component (FWHM = 1543 ± 34 km 
s _1 ) blueshifted by -157 ± 33 km s _1 with respect to Nl. This 
model is interpreted as unresolved rotation 



the narrow component in the rest frame. As in the previous 
case, it is possible that the observed narrow line splitting 
represents unresolved rotation, or traces a large-scale bipolar 
outflow in the gas. 

The very broad component is only detected in high 
ionization emission lines ([OIII]AA4959,5007, HeIIA4686 and 
NeVA3426). The other observed emission lines, with the ex- 
ception of H/3 and the Hcv+ [Nil] AA6549, 6583 complex, are 
adequately modeled using a model including the two nar- 
row components plus the broad component. In the case of 
the H/3 line and the Hq+[NII] complex the profiles the blue 
and red wings are not adequately fitted with a model com- 
prising the two narrow plus a broad component (see Figure 
I13p . Accounting for a fourth, blueshifted component does 
not improve the quality of the fit. In fact, when a forth, 
broad component consistent with that detected in [OIII] is 
included during the modelling, DIPSO returns negative flux 
values for this component. Therefore, we used the approach 
of including a fourth, unconstrained component (i.e. shift, 
width and flux of this component are free parameters) dur- 
ing the modelling of the H/3 and Ha+[NII] emission lines, 
and the results are shown in Figure [131 This fourth compo- 
nent has a FWHM of 2345± 368 km s' 1 and is redshifted by 
188 ± 100 km s _1 with respect to the narrow component at 
rest frame; its Hcv/H/3 ratio of 18±2.7 implies a redd ening of 
E(B-V)=1.68±0.14 assuming a lCalzetti etafl (120001 ) extinc- 
tion law. At the same time, the nuclear spectrum of the ob- 
ject shows an unusually red continuum compared with other 
local ULIRGs (Rodriguez Zaurin et al. 2009, see also Hill, G. 
et al. 1987). Overall, the properties of the nuclear spectrum 
of this source are consistent with it containing a moderately 
reddened Seyfert 1 nucleus. This is consistent with the cl ear 
detection at near-IR wavelengths bv lVeilleux et~afl (|l997l l of 
a broad Paa line of similar width (FWHM~1954 km s _1 ) 
to the broad Balmer lines detected in our optical spectrum. 

Note that a telluric absorption feature might affect the 
blue wing of the Ha+[NII] complex, as shown in Figure [141 
Nonetheless, we have successfully fitted the Hq+[NII] blend 



IRAS F23060 + 0505, Ha 

1 . 4 I 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r- 



1.2 - 

r i -° ~ 

S 0.8 - 

^ 0.6 - 

0.4 - 

1 0.2 - 

- 




6440 6480 6520 6560 6600 6640 6680 
Wavelength (A) 

Figure 14. Figure using the Ha+[NII] complex to illustrate 
the effects of the telluric absorption feature in the case of 
F23060+0505. Although the telluric feature affects the shape of 
the profile for the complex, it is still possible to model the He* and 
the [NII]AA6549,6583 emission lines using the [OIII] model. The 
red, solid line is our best fitting model. 



using a BLR component that has the same width and shift 
as that fitted to the H/3 line. 

Regarding the ionization mechanisms, we find that the 
line ratios measured for the three components observed in 
all the emission lines are consistent with Sy2 ionization. 
As in some previous cases, when only the shifts between 
the different components and their widths are fixed, the 
[SII]6716/6731 line ratio associated with N2 was higher than 
1.42 (the low density limit), while this value was lower 
than 0.44 (the high density limit) for the broad compo- 
nent. Following the approach described before, we set the 
values of this line ratio to be 0.44 and then change the value 
for each component individually in steps on 0.02 until the 
fit was deemed not adequate. Acceptable fits are obtained 
for any line ratio value within the allowed range (0.44 < 
[SII] 6716/6731 < 1.42). Finally, we note that the degener- 
acy issue associated with the [SII] AA6716, 6731 doublet does 
not change the results related to the ionization mechanism, 
i.e. Sy2-like line ratios are always found for all the different 
components. 

F23233+2817. The [OIII] model for this very compact ob- 
ject comprises 4 components: a narrow component (FWHM 
= 179 ± 4 km s" 1 ), an intermediate component (FWHM 
= 640 ± 24 km s _1 ) blueshifted by -161 ± 15 km s _1 
with respect to the narrow component, a broad component 
(FWHM = 1511 ± 118) blueshifted -447 ± 41 km s" 1 with 
respect to the narrow component, and a very broad compo- 
nent (FWHM = 2706 ± 422 km s^ 1 ) that is blueshifted by 
-3433 ± 904 km s _1 with respect to the narrow component. 
As shown in Figure [T] the flux contribution of the very broad 
component is relatively small. Therefore, the uncertainties 
associated with the corresponding shift and width are large. 
As in the previous object, this component is only required to 
obtain an adequate fit in the case of the [OIII] AA4959, 5007. 
All the other observed lines can be modelled with the first 
three components listed before. 

As in most of the cases, when only the shifts between 
the different components and their widths of are fixed, the 
[SII]6716/6731 line ratio associated with the intermediate 
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and broad components are outside the allowed range. There- 
fore, following the approach described above we set the ratio 
to 0.44 for the intermediate and broad components and then 
increase the value of each of these two components in steps 
of 0.02 until the fit was considered poor. In this way, we 
estimate an upper limit for the [SII]6716/6731 ratio of 1.0 
for these two components. 

Regarding the ionizing mechanisms, we find that the 
line ratios measured for the narrow component are consis- 
tent with NLAGN ionization, or fall close to the HII/Sy2 
limiting region in the diagrams. On the other hand, the line 
ratios measured for both the intermediate and the broad 
components are consistent with Sy2 ionization. Note that 
that a Sy2 spectral type is found for these two components 
independently of the [SII]6716/6731 line ratio used during 
the modelling. 

F23327+2913. Yuan et al. (2010) recently classified this 
double nucleus system as Sy2:, although it was originally 
classified as LINER bv lVeilleux et alj \ 19991) . For the work 
presented here we concentrate on t he 5kp c-I aperture de- 
scribed in IRodrfguez Zaurm et al.l (|2009l ) . This aperture 
covers the emission from the southern source, and is the 
only aperture for which the corresponding extracted spec- 
trum shows evidence for outflows. 

The main problem when trying to model the emission 
lines in this object is that the [OIII]AA4959,5007and H/3 
emission lines have low equivalent widths, and are poten- 
tially affected by noise and residual structure in the underly- 
ing continuum. Therefore these lines are not suitable for de- 
riving a model to fit to the other lines in the spectrum. As in 
the case of F00188-0856, we use the approach of first fitting 
the Ha+[NII] complex and then model the other emission 
lines observed in the spectrum with the same components 
and the same velocity widths and shifts as Ha. The results 
of the Ha+[NII] modelling are shown in Figure [T] 

The best fitting model comprises 2 components: a nar- 
row component (FWHM = 109 ± 2 km s _1 ) and an interme- 
diate component (FWHM = 805 ± 13 km s" 1 ) blueshifted 
by -92 ± 4 km s _1 with respect to the narrow component. 
The narrow component makes a negligible contribution to 
the flux of the H/3 emission line. In fact, the H/3 emission 
line can be modeled using only the intermediate component. 
Therefore, it is not possible to study the ionization mech- 
anism associated with that component. Regarding to the 
intermediate component, we find that the line ratios mea- 
sured for that component are consistent with NLAGN or 
HII ionization. 

F23389+0300. Two 5 kpc apertures (5kpc-I and 5kpc- 
II) were used during the investigation of the stellar 
populations of this do u ble n u cleus system presented in 
IRodrfguez Zaurm et all l|2009l , |2010T) However, aperture 
used for this study is the 5kpc-I aperture, centered on 
the northern sourc e, which includes the AGN nucleus (e.g. 
iNagar et all . 120031 ) 

The [OI] AA6300,6393 emission lines, which are stronger 
powerful than the observed [OIII] AA4959,5007 emission 
lines, are particularly impressive in this source. In terms 
of the modelling, it was not possible to find one model that 
adequately fits all the detected emission lines. In this partic- 
ular case, each emission line was fitted using an individual 
model. For completeness, Figure [T] shows the modelling re- 



sult for the [OIII] emission lines. The model shown in the 
figure comprises a narrow component (FWHM = 289 ± 7 
km s _1 ) plus a very broad component (FWHM = 2223 ± 30 
km s _1 ) that is redshifted by 47 ± 13 km s _1 with respect 
to the narrow component. 

As in the cases of F13451+1232, F14394+5332 and 
F16156+0146, the fact that we do not find one adequate 
model for all emission lines suggest that, instead of a discrete 
number of kinematic components, there may be a continu- 
ous range of density, ionization and kinematics through the 
Narrow Line Region (NLR), leading to each line having a 
different profile. No attempt has been made to estimate line 
ratios of the different kinematic components for this galaxy. 



